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A polydisperse system contains particles with at least one attribute a (such as particle size in 
colloids or chain length in polymers) which takes values in a continuous range. It therefore has an 
infinite number of conserved densities, described by a density distribution p{cr). The free energy 
depends on all details of p((t), making the analysis of phase equilibria in such systems intractable. 
However, in many (especially mean-field) models the excess free energy only depends on a finite 
number of (generalized) moments of p{cy); we call these models truncatable. We show, for these 
models, how to derive approximate expressions for the total free energy which only depend on such 
moment densities. Our treatment unifies and explores in detail two recent separate proposals by the 
authors for the construction of such moment free energies. We show that even though the moment 
free energy only depends on a finite number of density variables, it gives the same spinodals and 
critical points as the original free energy and also correctly locates the onset of phase coexistence. 
Results from the moment free energy for the coexistence of two or more phases occupying compa- 
rable volumes are only approximate, but can be refined arbitrarily by retaining additional moment 
densities. Applications to Flory-Huggins theory for length-polydisperse homopolymers, and for 
chemically polydisperse copolymers, show that the moment free energy approach is computationally 
robust and gives new geometrical insights into the thermodynamics of polydispersity. 
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I. INTRODUCTION 

The thermodynamics of mixtures of several chemical species is, since Gibbs, a well established subject (see e.g. [Q). 
But many systems arising in nature and in industry contain, for practical purposes, an infinite number of distinct, 
though similar, chemical species. Often these can be classified by a parameter, cr, say, which could be the chain length 
in a polymeric system, or the particle size in a colloid; both are routinely treated as continuous variables. In other 
cases (see e.g. |^-|^) c is instead a parameter distinguishing species of continuously varying chemical properties. In 
the most general case, several attributes may be required to distinguish the various particle species in the system 
(such as length and chemical composition in length-polydisperse random copolymers) and a is then a collection of 
parameters B . The thermodynamics of polydisperse systems (as defined above) is of crucial interest to wide areas of 
science and technology; it is sometimes also referred to as "continuous thermodynamics" (see e.g. |^). 

Standard thermodynamic procedures [|l| for constructing phase equilibria in a system of volume V containing M 
different species can be understood geometrically in terms of a free energy surface f{pj) (with / = F/V, and F 
the Helmoltz free energy) in the Af-dimensional space of number densities pj. Tangent planes to / define regions of 
coexistence, within which the free energy of the system is lowered by phase separation. The volumes of coexisting 
■ phases follow from the well-known "lever rule" [Q . Here "surface" and "plane" are used loosely, to denote manifolds 
of appropriate dimension. This procedure becomes unmanageable, both conceptually and numerically, in the limit 
Af — > oo which formally defines a polydisperse system. There is now a separate conserved density p{'j) for each value 
of cr; p((t) is in fact a density distribution and the overall number density of particles is written as p = Jdcr p{(t). The 
free energy surface is / = f[p{(j)] (a functional) which resides in an infinite-dimensional space. Gibbs' rule allows the 
coexistence of arbitrarily many thermodynamic phases. 

To make phase equilibria in polydisperse systems more accessible to both computation and physical intuition, it 
is clearly helpful to reduce the dimensionality of the problem. Theoretical work in this area has made significant 
progress in finding simplified forms of the conditions for phase coexistence, thus making these more numerically 
tractable | |:0-|l7|. Our aim is to achieve a similar reduction in dimensionality on a higher level, that of the free 
energy itself. We show that, for a large class of models, it is possible to construct a reduced free energy depending 
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on only a small number of variables. From this, meaningful information on phase equilibria can be extracted by 
the usual tangent plane procedure, with obvious benefits for a more intuitive understanding of the phase behavior 
of polydisperse systems. As an important side effect, this procedure also leads to robust algorithms for calculating 
polydisperse phase equilibria numerically. In particular, such algorithms can handle coexistence of more than two 
phases with relative ease compared to those used previously [[7| p^|JT5|Jl^ , [l8| -p5| . 

A clue to the choice of independent variables for such a reduced free energy comes from the representation of 
experimental data. We recall first the definition of a "cloud point" (see e.g. This is the point at which, 

for a system with a given density distribution /o(cr), phase separation first occurs as the temperature T or another 
external control parameter is varied. The corresponding incipient phase is called the "shadow" . Now consider diluting 
or concentrating the system, i.e., varying its overall density p while maintaining a fixed "shape" of polydispersity 
n{a) = p{a)/p. We will find it useful later to refer to the collection of all systems p{a) — pn{a) which can be obtained 
by this process as a "dilution line" (in the space of all density distributions p{a)). Plotting the cloud point temperature 
T versus density p then defines the cloud point curve (CPC), while plotting T versus the density of the shadow gives 
the shadow curve The differences in the shapes n{a) of the density distributions in the different phases are hidden 
in this representation; only the overall densities p in each phase are tracked. The density p is a particular moment 
of p(cr) (the zeroth one); higher order moments would be given by Jda a"^ p{a). Generalizing slightly, this suggests 
that our reduced free energy should depend on several (generalized) moment densities pi — Jda Wi{a)p{a), defined by 
appropriate (linearly independent) weight functions 'Wi{a). Ordinary power-law moment densities are included as the 
special case Wi{(j) — a^. Whatever the choice of moments, we insist on wq = 1 so that in all cases the zeroth moment 
density po coincides with the overall number density p. 

As an illustration, consider the simplest imaginable case where the true free energy / already has the required form, 
i.e., it depends only on a finite set of K moment densities of the density distribution: 

/-/(p,), t = l...K OT t^O...K~l (1) 

where the range of i depends on whether po = p is among the K moment densities on which / depends. (From now 
on, we use the notation i = 1 . . . K inclusively, to cover both possibilities.) In coexisting phases one demands equality 
of particle chemical potentials, defined as p{a) = 6f /6p{a), for all a. Because 

' ^ ^ i ' i 

this implies that all "moment" chemical potentials, p,i = df /dpi, are likewise equal among phases. The second 
requirement for phase coexistence is that the pressures or osmotic pressures [ p8| , 11, of all phases must also be equal. 
But from 



-11 = f - J dap,{a)p{a) = / - Y^PtPi 



one sees that this again involves only the moment densities pi and their chemical potentials fii. Finally, the moment 
densities also obey the "lever rule", as follows. Let the overall density distribution of a system of volume V be 
p^'^^(cr); we call this the "parent" distribution. If (after a lowering of temperature, for example) this parent has 
split into p coexisting "daughter" phases with a-distributions p^"^a), each occupying a fraction u*^"-' of the total 
volume (a = l...p), then particle conservation implies the usual lever rule (or material balance) among species: 
^^t;(")p(")(cr) ^ p(o)(cr),VCT. Multiplying this by a weight function Wi (a) and integrating over a shows that the lever 
rule also holds for the moment densities: 

f (2) 

a=l 

These results express the fact that any linear combination of conserved densities (a generalized moment density) 
is itself a conserved density in thermodynamics. We have shown, therefore, that if the free energy of the system 
depends only on K moment densities pi . . . pK we can view these as the densities of K "quasi-species" of particles, 
and construct the phase diagram via the usual construction of tangencies and the lever rule. Formally this has reduced 
the problem to finite dimensionality, although this is trivial here because /, by assumption, has no dependence on 
any variables other than the pi {i = 1 . . . K). 

Of course, it is uncommon for the free energy / to obey (|^). In particular, the entropy of an ideal mixture (or, 
for polymers, the Flory-Huggins entropy term) is definitely not of this form. On the other hand, in very many 
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thermodynamic (especially mean field) models the excess (i.e., non- ideal) part of the free energy does have the simple 
form (|^). In other words, if we decompose the free energy as (setting fee = 1) 

/ = -Tsid + f, Sid = - J da p{a) [Inp(cr) - 1] (3) 

then the excess free energy / is a function of some moment densities only: 

f = fip^), l^l...K (4) 

Examples of models of this kind include polydisperse hard spheres (within the generalization by Salacuse and Stell j|] 
of the BMCSL equation of state 1^^^ ), polydisperse homo- and copolymers and van der Waals fluids with 



factorized interaction parameters |10|] . With the exception of a brief discussion in Sec. Vl, this paper concerns models 
with free energies of the form (raM, which we call "truncatable" . (This terminology emphasizes that the number K 
of moment densities appearing in the excess free energy of truncatable models is finite, while for a non-truncatable 
model the excess free energy depends on all details of p(cr), corresponding to an infinite number of moment densities.) 
In what follows, we regard each model free energy as given, and do not discuss the issue of how good a description of 
the real system it offers, nor how or whether it can be derived from an underlying microscopic Hamiltonian. Whenever 
we refer to "exact" results, we mean the exact thermodynamics of such a model as specified by its free energy. 
Different authors refer in different terms to Sid in (^. Throughout this paper will call the quantity 

Sid = - / da p{a) [lnp{a) - 1] 



(which is, up to a factor of — T, the ideal part of the free energy density) the "entropy of an ideal mixture" . By 
writing p{a) — pn{a), where n{a) is the normalized distribution of a, this can be decomposed as 



Sid 



—p{\np — l) — p J da n{a) \n n{a) (5) 



The first term on the right hand side is the "entropy of an ideal gas", while the second term gives the "entropy of 
mixing" . The prefactor p reflects the fact that both are expressed per unit volume rather than per particle. 

In principle, the entire phase equilibria for any truncatable system (obeying (^) can be computed exactly by a 
finite algorithm. Specifically, the spinodal stability criterion involves a X-dimensional square matrix [pTl- lil w hereas 
calculation of p-phase equilibrium involves solution of {p— 1)(-R' + 1) coupled nonlinear equations (see Sec. HIE). This 
method has certainly proved useful [p|-^|j9| JlO| , p^ , but is cumbersome, particularly if one is interested mainly in cloud- 
and shadow-curves, rather than coexisting compositions deep within multi-phase regions ||^,^,^, 26|. Various ways of 
simplifying the procedure exist [^11 15|-p^, but there has previously not been a systematic alternative to the full 
computation. Note also that the nonlinear phase equilibrium equations permit no simple geometrical interpretation 
or qualitative insight akin to the familiar rules for constructing phase diagrams from the free energy surface of a finite 
mixture. This further motivates the search for a description of polydisperse phase equilibria in terms of a reduced 
free energy which depends on only a small number of density variables. 

In previous work, the authors originally arrived independently at two definitions of a reduced free energy in terms 
of moment densities (moment free energy for short) |31p^ ]. Though based on distinctly different principles, the two 
approaches led to very similar results. We explain this somewhat surprising fact in the present work; at the same 
time, we describe the two methods in more detail and explore the relationship between them. We also discuss issues 
related to practical applications and give a number of example results for simple model systems. 

The first route to a moment free energy takes as its starting point the conventional form of the ideal part of the 
free energy, /id — T ^ da p{a)\h\ p{a) — 1], as in (^). This can be thought as defining a (hyper-) surfa ce, with the 



"horizontal" co-ordinate being p{a) and the "height" of the surface /id[jo((T)]. As explained in Sec. II A , this surface 
can then be projected geometrically onto one with only a finite set of horizontal co-ordinates: these are chosen to be 
the moment densities appearing in the excess free energy. Physically, this corresponds to minimizing the free energy 
with respect to all degrees of freedom in p{a) that do not affect these moment densities. We call this first route the 
"projection" method. 

The second approach rederives the entropy of mixing in the ideal part of the free energy in a form that depends 



explicitly only on the chosen moment densities. As described in Sec. II B , the expression that results is intractable 
in general because it still contains the full complexity of the problem. However, in situations where there are only 
infinitesimal amounts of all but one of the phases in the systems, the entropy of mixing can be evaluated in a closed 
form. Applying this functional form in regimes where it is not strictly valid (i.e., when phases of comparable volumes 
coexist) generates a moment free energy by this "combinatorial" method. 
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We show the equivalence of the two approaches in Sec. ilC. There, we first demonstrate that the general form 
of the entropy of mixing obtained by the combinatorial method can be transformed to the standard expression 
— Jdan{a) hin{a). Second, we show that the moment free energies arrived at by the two methods are in fact equal, 
with the pro jection method being slightly more generally applicable. 

In Sec. Ill, we then discuss the properties of the moment free energy (as obtained from either the projection or the 
combinatorial method). By construction, it depends only on a (finite) number of moment densities, and so achieves the 
desired reduction in dimensionality. The moment densities can be treated as densities of "quasi-species" of particles, 
and the standard procedures of the thermodynamics of finite mixtures can be applied to the moment free energy 
to calculate phase equilibria. This is only useful, however, if the results are faithful to the actual phase behavior 
of the underlying polydisperse system (as modeled by the given truncatable free energy). We show that this is so: 
In fact, the construction of our moment free energy is such that exact binodals (cloud-point and shadow curves), 
critical (and multi-critical) points and spinodals are obtained. Beyond the onset of phase separation, where coexisting 
phases occupy comparable volumes, the results are not exact, but can be refined arbitrarily by adding extra moment 
densities. This procedure is necessary also to ensure that, in regions of multi-phase coexistence, the correct number of 
phases is found. In Sec. we discuss the practical implementation of our method, followed by a number of examples 
(Sec. Our results are summarized in Sec. where we also outline perspectives for future work. 

Note that throughout this work, we focus on the case of phase coexistence at fixed volume. However, as described in 
App. the formalism can easily be applied to scenarios where the (mechanical or osmotic) pressure is fixed instead. 



II. DERIVATION OF MOMENT FREE ENERGY 



A. Projection method 

The starting point for this method is the decomposition (^) for the free energy of truncatable systems (^): 



f = T jdapia) 



In- 



fiP^) (6) 



As explained in the previous section, truncatable here means that the excess free energy / depends only on K moment 
densities pi. Note that, in the first (ideal) term of (^, we have included a dimensional factor R{a) inside the logarithm. 
This is equivalent to subtracting T Jda p{a) In R{a) from the free energy. Since this term is linear in densities, it has 
no effect on the exact thermodynamics; it contributes harmless additive constants to the chemical potentials /^(cr). 
However, in the projection route to a moment free energy, it will play a central role. 

We now argue that the most important moment densities to treat correctly are those that actually appear in 
the excess free energy f{pi). Accordingly we divide the infinite-dimensional space of density distributions into two 
complementary subspaces: a "moment subspace" , which contains all the degrees of freedom of p{a) that contribute 
to the moment densities pi, and a "transverse subspace" which contains all remaining degrees of freedom (those that 
can be varied without affecting the chosen moment densities pi) Q. Physically, it is reasonable to expect that these 
"leftover" degrees of freedom play a subsidiary role in the phase equilibria of the system, a view justified a posteriori 
below. Accordingly, we now allow violations of the lever rule, so long as these occur solely in the transverse space. 
This means that the phase splits that we calculate using this approach obey the lever rule for the moment densities. 



but are allowed to violate it in other details of the density distribution p{a) (see Fig. ^in Sec. VB 1 for an example). 
These "transverse" degrees of freedom, instead of obeying the strict particle conservation laws, are chosen so as to 
minimize the free energy: they are treated as "annealed". If, as assumed above, / = fipi) only depends on the 
moment densities retained, this amounts to maximizing the ideal mixture entropy (sid) in (o), while holding fixed the 
values of the moment densities pi. Note that we are allowed, if we wish, to include among the retained densities 
"redundant" moments on which / has a null dependence. We will have occasion to do this later on. 

At this point, the factor R{a) in (|^), which is immaterial if all conservation laws are strictly obeyed, becomes 
central. Indeed, maximizing the entropy over all distributions p(a) at fixed moment densities pi yields 

(7) 



p(cr) = R{a) exp XiWi{a) 



where the Lagrange multipliers Xi are chosen to give the desired moment densities 

(8) 



= J da {a) R{a) eyip XiW,{a)^ 
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The corresponding minimum value of / then defines our projected free energy as a function of the moment densities 
Pi- 



i 

Here Spi- is the projected entropy of an ideal mixture. The first term appearing in it, po — Jda p{a), is the "zeroth 
moment" , which is identical to the overall particle density p defined previously. If this is among the moment densities 
used for the projection (or more generally, if it is a linear combination of them), then the term —TpQ is simply a linear 
contribution to the projected free energy /pr(pi) and can be dropped because it does not affect phase equilibrium 
calculations. Otherwise, po needs to be expressed - via the - as a function of the pi and its contribution cannot be 
ignored. We will see an example of this in Sec. ^ 

The projection method yields a free energy /pr(pO which only depends on the chosen set of moment densities. These 
can now be viewed as densities of 'quasi-species' of particles, and a finite-dimensional phase diagram can be constructed 
from /pr according to the usual rules, ignoring the underlying polydisperse nature of the system. Obviously, though, 
the results now depend on R{a) which is formally a "prior distribution" for the entropy maximization. To understand 
its role, we recall that the projected free energy fpripi) was constructed as the minimum of /[p(o')] at fixed pf, that 
is, /pr is the lower envelope of the projection of / onto the moment subspace. Crucially, the shape of this envelope 
depends on how, by choosing a particular prior distribution R{a), we "tilt" the infinite-dimensional free energy surface 
before projecting it. This geometrical point of view is illustrated, for a mixture of only two species, in Fig. |l|. 

To understand the effect of the prior R{cr) physically, we note that the projected free energy is simply the free 
energy of phases of a system in which the density distributions p{a) are of the form (0). The prior R{cr) determines 
which distributions lie within this "maximum entropy family" (or "family" for short) , and it is the properties of phases 
with these distributions that the projected free energy represents. Typically, one is interested in a system where a 
fixed overall "parent" (or "feed") distribution p^'^''{a) becomes subject to separation into various phases. In such 
circumstances, we should generally choose this parent distribution as our prior, R{a) = p^''^(cr), thereby guaranteeing 
that it is contained within the family (^. Having done this, we note that the projection procedure will be exactly 
valid, to whatever extent the density distributions actually arising in the various coexisting phases of the system under 
study are members of the corresponding family 



p(a)=pW(a)exp(5]A.w;,(a)) . (10) 



In fact, the condition just described holds whenever all but one of a set of coexisting phases are of infinitesimal volume 
compared to the majority phase. This is because the density distribution, p'^*'^(ct), of the majority phase is negligibly 
perturbed, whereas that in each minority phase differs from this by a Gibbs-Boltzmann factor, of exactly the form 
required for (p^); we show this formally in Sec. [II. Accordingly, our projection method yields exact cloud point and 



shadow curves. By the same argument, critical points (which in fact lie at the intersection of these two curves) are 
exactly determined; the same is true for tricritical and higher order critical points. Finally, spinodals are also found 



exactly. We defer explicit proofs of these statements to Sec. III. 

The projection method does, however, give only approximate results for coexistences involving finite amounts of 
different phases. This is because linear combinations of different density distributions from the family (|l^), corre- 
sponding to two (or more) phases arising from the same parent p*-"-' (a) , do not necessarily add to recover the parent 
distribution itself: 

Y^v^'^^p^'\a)e^p (^X^^w^a)^ ^ p^^Ha) (11) 

unless all except one of the w^") are infinitesimal. Moreover, according to Gibbs' phase rule, a projected free energy 
depending on n moment densities will not normally predict more than n+1 coexisting phases, whereas a polydisperse 



system can in principle separate into an arbitrary number of phases. We explain in Sec. Ill how both of these 
shortcomings can be overcome by systematically including extra moment densities within the projection procedure. 
How quickly convergence to the exact results occurs depends on the choice of weight functions for the extra moment 
densities; we discuss this point further in the context of the examples of Sec. 0. 
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FIG. 1. Illustration of free energy projection. The full infinite-dimensional free energy surface f[p{<y)] for a truly polydisperse 
system cannot be represented graphically, so a system with only two particle species is used instead. A single phase is then 
described by its density pair {ri,r2). (The projection method remains applicable, although the polydispersity parameter a is 
replaced by a discrete species index a £ {1, 2}, and integrals J da p{a) ... by Va ■ ■ ■) To illustrate the method we choose a 
model with an excess free energy depending on a single moment density pi = L\r\ + L2r2 (in fact Flory-Huggins theory for 
length polydispersity, discussed in more detail in Sec. see (|6^)). Both the transverse and the moment subspace are then 
one- dimensional, (a) 3d- view of free energy surface. The circle marks the position of the parent (rj"', rj"'); the thick line 
traces out the free energy of density pairs (ri,r2) in the corresponding family (^o|). (b) Top view, showing as dashed lines the 
transverse directions along which the moment density pi is constant, (c) Side view of the free energy surface, "looking down" 
the transverse direction, (d) The projected free energy is simply the lower envelope of the free energy surface as seen from this 
transverse direction, (e-h) The same free energy surface, now tilted differently by using a different parent (or prior). While this 
does not affect the location of any double tangent planes drawn to the full surface, it does produce a different "lower envelope" 
and hence a different projected free energy. For continuum polydispersity, the role of the density pair {r\,r2) is played by a 
density distribution p[a) of a continuous parameter a. 
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B. Combinatorial method 



Now we turn to the second method of constructing a moment free energy. As before, we recognize from the outset 
that the physics of the problem is contained in the excess free energy and therefore the most important variables 
are the moment densities which feature in this. On the other hand, the ideal part of the free energy (essentially the 
entropy of mixing, up to the trivial ideal gas term) is a book-keeping device that accounts for the number of ways 
of partitioning a priori a given size distribution between two or more coexisting phases. (To focus the discussion, 
we identify a as particle size in this section; the arguments do of course remain valid for a general polydispersity 
parameter a.) The entropy of mixing has its origin in factors in the partition function (where N is a, particle 

number) which are usually derived from the classical limit of quantum statistics. The origin of the factors in 

classical statistical mechanics has been the subject of debate ever since Gibbs . The problem is also connected with 
the question of the extensivity of entropy, and the Gibbs paradox. Below we present a completely classical derivation 
of the entropy of mixing which has the additional benefit of indicating the appropriate generalisation for moment 
densities. To start with, follow Gibbs and define a non-extensive free energy F' as a configuration space (T) integral 

e-^'/T^ f dVe-"'^. (12) 



No appears in this as all particles are distinguishable in classical statistical mechanics. 

Now consider two phases in coexistence as one joint system. Assume they occupy volumes V'^^^ and V^'^\ respec- 
tively, and contain N'^^^ and A^^^^ particles. Following literally the prescription in eq. (|l2|), the free energy of the joint 
system of A^^^^ -f A^^^^ = A^ particles is found from 



-F' /T _ (_F'(i)+F'(^')/T 



prtns 



(13) 



where the configuration space integral has been done in two parts. Firstly, for each way of partitioning the particles 
between the two phases, the individual configuration space integrals give a product of the individual partition functions. 
Secondly, and crucially, one must sum over the A^!/A^'^^^!A^'^^^! ways of partitioning the particles between phases. Now 
define an extensive free energy by reinserting the 1/A^! as though the particles were indistinguishable: 



e-^'"/^ = ^ /rfre-^/^, (14) 



Eq. ( [l3| ) can then be written as 



)prtns (15) 



where the average is taken over all partitions, with equal a priori probabilities. This result is the cornerstone of the 
combinatorial method. If we separate Fi„ into the entropy of an ideal gas and the remaining non-ideal (excess) parts, 
it is written as 

Fi„ ^ NT[\n{N/V) -1]+F (16) 

with F being the excess free energy. Note that F„^ does not yet contain an entropy of mixing term because it treats the 
particles as though they are indistinguishable. The entropy of mixing will reappear when we come to do the average 
over partitions in (|l^) . Note also that conventional thermodynamics follows from (|l^) if all particles are identical as 
far as their mutual interactions are concerned. There is then no entropy of mixing to consider, so that F^^ is just the 
conventional free energy F; and since the average over partitions is trivial one has F = + Similarly, as we 

show in Sec. [IC, one can recover from ( [l5| ) the conventional form of the entropy of mixing as given in (|^). 



We now show how the average over partitions in ( |l5|) results in an expression for the entropy of mixing which 
depends explicitly on moment variables. The key is to note that, by our assumption that we are dealing with a 
truncatable model, the excess free energy F = V f{pi) depends only on a limited number of moment densities pi (as 
well as volume V , temperature T, and other state variables which we suppress below). If the density po is among 
the moment densities pi - which we assume throughout this section - then specifying V and the pi is equivalent to 
specifying A^, V and the normalized moments rrii = pi/po {i > 0). So we can think of as a function of A^, V and 
the rrii. From (p^, Fm then depends on the same set of variables. The choice of the m^, which are moments of the 
normalized particle size distribution n{a), as independent variables is natural in the present context because we are 
considering the particle number A^ and volume V of the coexisting phases to be fixed. 
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To avoid notational complexity, and without loss of generality, we now specialize to the case where there is only 
one normalized moment m^, a generalized mean size 



m — J da w{a)n(a) 

We shall also write pm = pni — Jda w{a)p{a) for the corresponding moment density, and m^^^ and to'^' for the values 
of the generalized mean size in the first and second phase. The overall ("parent") size distribution is n'-^^{a) with 
mean to^''^ . In each of the two phases, in which and V are held fixed, Fin only depends on m, and so the average 
over partitions becomes 



Here V{m'--^">) is the probability distribution for the generalized mean size in the first phase, taken over partitions with 
fixed N'^^^ and iV'^^ with equal a priori probabilities. Note that given m'-^\ wS"^^ is fixed in the second phase by the 
moment equivalent of particle conservation: N^^^m'^-^^ + iV'^^m'^^^ = Nm^'^\ The integral in ( pT| ) can be replaced by 
the maximum of the integrand in the thermodynamic limit, because ln'P(m^^-') is an extensive quantity. Introducing a 
Lagrange multiplier for the above moment constraint then shows that the quantity pm has the same status as the 
density p = pQ itself: both are thermodynamic density variables. This reinforces the discussion in the introduction, 
where we showed that moment densities can be regarded as densities of "quasi-species" of particles. 

In what follows, we only need to refer to to'-'^' (not m^^-*) and therefore drop the superscript for brevity. Since 
In ■p(m) is extensive, we can write In'P(TO) = Ns{m) where s is the entropy of mixing per particle expressed as a 
function of the moment m (or, more generally, of the full set of normalized moments). The total free energy of the 
system then takes the form 

Fin - min f/^^ + F^f ^ - NTs{m), (18) 

m 

and we need to calculate the entropy of mixing s{m). We recognize that this quantity depends not only on the 
generalized mean size, but also on x = iV'-'^'/Af, the fraction of particles in the first phase. A formal result for s(m) 
can be obtained by first calculating the joint probability distribution V{xm,x) of xm and x. From this, we can find 
'P(m) according to |3|] 

^ V{^rn,x) ^ r{xm,x) 
^ ' Jdm'P{xm',x) V{x)/x ^ ' 

Writing the log probabilities in ( p!9| ) as lnP(a;m,a;) = Ns{xm,x) and hiVix) = Ns{x), we then have in the thermo- 
dynamic limit 

s{m) — s{xm, x) — s{x) (20) 
The method for calculating Vixm, x) proceeds by introducing an indicator function for each particle, deemed 



to be 1 if the particle is in the first phase, and if the particle is in the second |36 . Then Nx — 



Nxm — X]i3=i ^i'^i^i)- We now write the moment generating function for V{xm, m) as follows 



j d{xm) dx V{xm,x) exip[N{9x + Xxm)]— ^expj^fl^^-^ + A^^-^ eiw{<7i) 

N N 

'T\ exp{e,[e + Xw{<j,)]}) = 2-^ n{l + cxp[0 + Xw{a,)]} 

\ / prtns 

i=l ^ i=l 

= 2-^ exp ln{l + exp[e' + Xw{a^)]}j . 



(21) 



In the second equality we have used the fact that the e.; are independent discrete random variables taking the values 
and 1 with a priori equal probabilities. Taking logarithms of (|21) and dividing by N, we obtain on the right hand side 
an average over i — 1 ... of a function of ai. Since the ai are drawn independently from n^^'>{a), the law of large 
numbers guarantees that this average tends to the corresponding average over n^^^cr) in the limit N ^ oo (compare 
the discussion in Ref. |0). We therefore have the final result 
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-i-ln Jd{xm)dx exp{N[s{xm,x) + 9x + Xxm]} = J da n^°\a) ln{l + cxp[6 + Xw (a)]} (22) 

where we have dropped an unimportant additive constant (—In 2) and have replaced V{xm,x) by exp[Ns{xm,x)]. 
Now, the left hand side can be evaluated, by the steepest descent method, in the limit TV — > oo. Defining 

V'(6',A) jdan^^\a) ln{l + exp[6' + Au;(cr)]} 

we get s{xm, x) + Ox + Xxm ~ ijj{0, A) at the point where ds/dx + 6 — ds/d{xm) + A = 0. We recognize that the 
relationship between s{x,xm) and ip{O^X) is a double Legendre transform. Inverting this shows that 

s{xm, x) — ?A(0, X) ^ 9x ~ Xxm (23) 

where 

x^di^/de xm = dtp/dX (24) 

Eq. ( p3| ) is essentially the desired result; to obtain s(m) from (|20|), only s{x) remains to be determined. It is obtained 
from V{x) via exp[A'^s(x)] — V(x) — ^ d{xm)V{xm,x) which by inspection of the moment generating function is 
seen to correspond to the point A = 0. At this point, if) = ln(l + e^) and x = dij^/dO = e^/(l + e^). We find that 
= ln[x/(l — x)] and after some manipulation, 

s{x) — Tp ^ x9 ~ —a; In a; — (1 — a;) ln(l — x) (25) 

This is recognized as the standard entropy of mixing that is lost when a total of N particles is partitioned into 
A^(i) = xN and N^^^^ = (1 — x)N particles in two phases. The simple form of the result shows that the calculation 
leading to ([2^) is essentially a gen eralization of the Stirling approximation. 
Inserting (p3[) and (|25| ) into (|2C|), we finally have the desired result for s{m), 

s(m) — -iplO, X) ^ 9 X ~ X xm + a;lna; + (1 — a;) ln(l — a;) (26) 

Unfortunately this is mainly formal because neither the integral defining -0 nor the Legendre transform are likely to 



be tractable. However, we show in Sec. [IC that (|2^) is equivalent to the more conventional form of the entropy 
of mixing as given by the second term in in eq. (^. From a conceptual point of view, it should be noted that the 
conventional form is normally derived by "binning" the distribution of particle sizes a and taking the number of bins 
to infinity after the thermodynamic limit N ^ oo has been taken (see e.g. |^). In our above derivation, on the other 
hand, we have assumed that even for finite N all particles have different sizes (7^, drawn randomly from n^^\a)] the 
"polydisperse limit" is thus taken simultaneously with the thermodynamic limit. The relation between these two 
approaches - which lead to the same results - has been discussed in detail by Salacuse . Note that the first limit is 
physically more plausible for many homopolymer systems (where there may only be thousands or millions of species, 
with many particles of each) whereas the second limit is more natural for colloidal materials (and also some random 
copolymers) in which no two particles present are exactly alike, even in a sample of macroscopic size. 

Although the full result (p9) is intractable, progress can be made for a; — > 0, when the number of particles in one 
phase is much smaller than in the other. The limit a; — > implies 9 — » — oo and hence V'(^; ^) ^ exp[0 + /i(A)] 
where h — \iy Jdcm^"-' ((T)e^'^^°'-' is a generalized cumulant generating function for 'n!'^\a). From this we derive 
X = d'ip/d9 — -0, giving 9 — Inx — h. The generalized mean size is given by xm = dijj/dX = ipdh/dX and hence 
m = dh/dX. Inserting these results into ( p3| ) then shows that 

s(a;m,a;) = a; — a;lna; + a;(ft. — Am), ttt 

oX 

and hence from ( |2^ ) 

s{m)=x{h-Xm) + Oix'^). (27) 

Noting that a;A^ = N^^'' is the number of particles in the small phase, the total free energy (|l^ ) can then be written 
as 

i^in = F^;^^ + - N^'^T{h - Xm) (28) 
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The term {h — Am), being multipUed by — iV^^^T, can now be interpreted as the entropy of mixing per particle in 
the small phase. It arises from the deviation of the (generahzed) mean size m = m^^^ in the smaU phase from the 
mean size m^^^ of the parent, and is given by the Legendre transform of the generahzed cumulant generating function: 

s^i^(m) = /i - Am, where "i = ^ = ^^^ j dan'^°\a) e^'^'-'^l (29) 

Let us now examine the relation between Smix(™), thereby defined, and s(m) introduced previously. By construction, 
Ns{m) is the entropy of mixing of the system as a whole; both phases contribute to this. The result (^^, which can 
also be written as 

Ns{m) = 7V(i)s,„i^(m) (30) 

appears to contradict this, seeing as it does not contain a term proportional to TV^^^ . The resolution of this paradox 
comes from the neglected 0{x^) terms in (^7|). In fact, using the single-phase entropy of mixing defined in eq. (^9|) - 
and reinstating on the r.h.s. the superscript on m = m*^^^ - we can write ( pO| ) in the more symmetric form 

Ns{m) = A^(i)s„,ix(m(i)) +iV(2)s„,i,(m(2)) 

This is still correct to leading (linear) order in x because the added term is 0{x'^). [To see this, use the fact that 
m^^-* = m^^"^ — xm^^''/{l — x) and hence m^^^ — m*^"^ = 0{x). From ([29|) it can then be deduced that Smix(m'^^-') ~ 



(m^ ' — ') = 0{x ).] Similarly, we can rewrite the total free energy (28) as 



Fi„ - [4^) - iV(i)Ts„,ix(m(2))] + [4^) _ iv(2)j^,^^.^(^(2))] (31) 

This shows that in the limit where one of the two phases is much larger than the other one, the free energy of each 
of the phases - now including the entropy of mixing - is given by 

F - Fin - NTs,,,i^{m) ^ -NT [- HN/V) + 1 + s„ix(m)] + F (32) 

This expression, which only depends on the particle size distribution through the generalized mean size m, is our 
desired moment free energy. 

The Legendre transform result (|2g|) for the (single-phase) entropy of mixing is appealingly simple; we will illustrate 



its application in Sec. VA, using polydisperse Flory-Huggins theory as an example. As discussed briefly in App. 
eq. (^9|) also establishes an interesting connection to large deviation theory. However, the most important aspect 
of the result ( p9| , ^2|) is that - as we have shown ~ it gives exact results for the limiting case N^^^^ <^ iV^^^^ (that is, 
X ^ 0). This includes two important classes of problems, that are also handled exactly by the projection method 
(for essentially the same reason). The first is the determination of spinodal curves and critical points. Intuitively 
these are exact because they are related to the stability of the system with respect to small variations in its density 
or composition, which can be probed by allowing fluctuations to take place in a vanishingly small subregion. The 
second concerns the cloud point and shadow curves. Again intuitively, these are exact because by definition only an 



infinitesimal amount of a second phase has appeared. Formal proofs of these statements are given in Sec. [II. 

If the result (^,^,^ is applied in the regime where it is no longer strictly valid, i.e., where the two coexisting 
phases contain comparable number of particles (or, equivalently, occupy comparable volumes), then approximate 
two-phase coexistences can be calculated. It is also straightforward to show that the analog of eqs. (|29|j33,||) holds 
in the case where two or more small phases coexist with a much larger phase; application in the regime of comparable 
phase volumes then provides approximate results for multi-phase coexistence. 



To end this section, let us state in full the analogue of (E9 32) for the case of several moment densities, restoring 



the notation used in the previous sections. The square bracket on the r.h.s. of (BSf) is the moment expression for the 



entropy of an ideal mixture. If, as in Sec. 11 A, we measure this entropy per unit volume (rather than per particle, as 
previously in the current Section) and generalize to several moment densities, we find by the combinatorial approach 
the following moment free energy: 

/comb = -T'Scomb + /, Scomb = -Po(lnPO - + PO \ h - XiTTli , (33) 



with 



h — hi I d(7 71^°^ (cr) exp ^ XiWi {a) 



(34) 



Here po = N/V is the particle density as before. 
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C. Relation between the two methods 



Although eq. ( p3| ) looks somewhat different from the corresponding result (||) for /pr obtained by the projection 
method, the two methods are, mathematically, almost equivalent, as we now show. (For brevity we return to the 
case of a single moment m and continue to refer to the polydisperse feature a as 'size'.) First, consider the exact 
expression for the entropy of mixing (|2^) derived within the combinatorial method. This should be equivalent to the 
conventional result used as the starting point (|^) in the projection method. To see this, write the Legendre transform 



conditions (24) out explicitly: 



1 + exp[6' + Aw(cr)] ' J ^ ^ ^ ^ 1 + exp[6' + Aw(cr)] 

Here we have reinstated the superscript to show that rrS^^ is the value in phase one of the (generalized) moment m 

to(i) = / da w{a)n'' ^\a) 



Comparing with ( |35| ) , we can identify the (normalized) particle size distribution in that phase as ||38 

(1) nW(a) exp[g + Au;(a)] 
'(a) — 



X 1 + exp[0 + Aw(cr)] 

Particle conservation xrv- i)(cr) + (l-a;)n(2)(cr) {a) implies a similar form for the size distribution in phase two: 



I- X 1 + exp[6' + \w{a)] ' 

It then takes only a few lines of algebra to show that the entropy of mixing ( p6| ) can be written (for any x) as 

Sim) ^-xjda (a) In ^ - [1 - x) J da n^^ (a) In ^ (36) 
^-x fdan^^^ (a) In n^^^ (a) - (1 - x) [da n^^) (a) In n^^) (a) + f da n^"^ (a) In n^o) (a) (37) 



The second equation, with the third (constant) term discarded, is the standard result (compare eqs. (^J^ jAl|) ). However 
the first expression appears more natural, and is retained, within the combinatorial derivation; it embodies the 
intuitively reasonable prescription that entropy is best measured relative to the parent distribution n'^^\a). This 



was also an essential ingredient of the projection method, as described in Sec. HA, where the "prior" in the entropy 
expression was likewise identified with the parent. Retaining the parent as prior also avoids subtleties with the 
definition of the integrals in (^^ in the case where the phases contain monodisperse components, corresponding to 
(5-peaks in the density distributions 

Having established that the rigorous starting points of the combinatorial and projection method are closely related, 
we now show that the subsequent approximations also lead to essentially the same results. The relevant approximations 
are use of ( ^3|) for /comb (which adopts the small x form for Smix) in the combinatorial case and use of (|) for /p, 
(which minimizes over transverse degrees of freedom) in the projection case. First note that the density po of phases 



within the family (10) is given by 

Pa = J da p'^^\a) exp XiWi{a)j 

Comparing this with (34), and using the fact that n^'^\a) = p^^\a)/ p'"q \ one sees that 

PO = y dap(")(a)exp |^ A,;u;,(a) ) = p<,"^e^°+" 



where the sum over i now includes i = 0. Solving for h, one has h — — Ao + \n[pQ/ p^^^). The ideal mixture entropy 
derived by the combinatorial route, eq. (^), can thus be rewritten as 
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This is identical to the projected entropy Spr, eq. except for the last term. But by construction, the combinatorial 
entropy assumes that po ~ the overall density - is among the moment densities retained in the moment free energy. 
The difference Scomb — Spr = ^Po In Po"'' is then linear in this density, and the combinatorial and projection methods 
therefore predict exactly the same phase behavior. 

In summary, we have shown that the projection and combinatorial methods for obtaining moment free energies 
give equivalent results. The only difference between the two approaches is that within the projection approach, one 
need not necessarily retain the zeroth moment, which is the overall density po = a-s one of the moment densities on 
which the moment free energy depends. If po does not appear in the excess free energy, this reduces the minimum 
number of independent variables of the moment free energy by one (see Sec. M for an example). 



III. PROPERTIES OF THE MOMENT FREE ENERGY 



In the previous sections, we have derived by two different routes (namely (|9|) and (p3[)) our moment free energy 
for truncatable polydisperse systems. We now investigate the properties of this moment free energy, which we 
henceforward denote /m and write in the form (|^) 

fmiPi) = -TSaiipi) + f{Pi), Sm = PO - ^ >^iPi ] (38) 

i 

here Sm is the "moment entropy" of an ideal mixture. In particular, we formally compare the phase behavior predicted 
from this free energy (by treating the pi as densities of "quasi-species" of particles, and applying the usual tangency 
construction) with that obtained from the exact free energy of the underlying (truncatable) model 

f[p{<j)] I dap{a) [\np{a) - 1] + /(p,). (39) 



Let us first collect a few simple properties of the moment free energy (^8|) which will be useful later. Recall that ( p3q ) 
faithfully represents the free energy density of any phase with density distribution in the family 



p(a)=p(o)(a)exp|^^A,u;,(a)^ 



(40) 



where p^*'-'(cr) is the density distribution of the parent. The moment densities pi are then related to the Lagrange 
multipliers \i by 

p^ = Jd<Tw^{cr) p^°\a)exp ^^XjWj{(T)j (41) 

If we regard po as a function of the A^, then from ( plj ) we have dpo/dXi = pi. Together with eq. (^8|), this implies that 
the moment entropy Sm has the structure of a Legendre transform p9| . For the first derivatives of Sm with respect to 
the moment densities, this yields 

dsn^ _ _^ 
dpt 

while the matrix of second derivatives is the negative inverse of the matrix of "second-order moment densities" pij [|oj : 
^^ = -(M-i),„ (M),, = ^^= y'da^.(a)«;,(a)p(°)(a)exp|^^Afeu;,(a)j ^p,,. (42) 

The chemical potentials pi conjugate to the moment densities follow as 

p^ = ^=TX,, + ^ = TX, + p, (43) 
opi api 
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and their derivatives w.r.t. the pi, which give the curvature of the moment free energy, are 

TT - = nM-% + (44) 

dpj opidpj opiopj 

The pressure, finally, is given by 



Eqs. ( |43[445| ) are all calculated via the moment free energy. The three corresponding quantities obtained from the 
exact free energy ( |3^ ) are, first, the chemical potentials conjugate to p{(j): 

Pi.^) = =T^'^Pi^) (46) 

5p{a) ^ 

second, their derivatives w.r.t. p(cr): 

^ix(a) S^f TSja -a') , ^ d^f 

,s ^ ..... ,^ = -r~\ ^/ a a W^(cr)Wj(a) (47) 

and, third, the resulting expression for the pressure: 

'U = f-Jda /i(a)p(a) = -Tpo + f 'Y. P^P^' (48) 

i 

The last of these is identical to the result (|4|) derived from the moment free energy. 

We note one important consequence of the form of the exact chemical potentials (^6|) for truncatable models: If 
two phases p'^^^a) and p^^^(tT) have the same chemical potentials p{cr), then the ratio of their density distributions 
can be written as a Gibbs-Boltzmann factor: 



pW(a) _ 

p(2)(cT) " P 



(49) 



This implies that if one of the density distributions is in the family ([40|), then so is the other. The same argument 
obviously applies if there are several phases with equal chemical potentials. Conversely, we have for the chemical 
potential difference between any two phases in the family ( ^0| ) 

Ap{<j) =J2iT^K + Ap,)w,{<j) ^J2Ap,w,{a). (50) 

i i 

The chemical potentials p{(j) of two such phases are therefore equal if and only if their moment chemical potentials 
Pi are equal. Combining this with the fact that the pressure (|4g) derived from the moment free energy is exact 
(compare (|4^)), we conclude: Any set of (two or more) coexisting phases calculated from the moment free energy 
obeys the exact phase equilibrium conditions. That is, if they were brought into contact, they would genuinely coexist. 
(However, they will not necessarily obey the lever rule with respect to the parent p^^^ (a) of the given family, unless 
all but one of them - which then coincides with the parent - are of vanishing volume; see (111]).) 



A. General criteria for spinodals and (multi-) critical points 

We now demonstrate that, for any truncatable model, the moment free energy gives exact spinodals and (multi-) 
critical points. By this we mean that, using the moment free energy, the values of external control parameters - such 
as temperature - at which the parent phase p^^^a) becomes unstable or critical can be exactly determined. Our 
argument treats spinodals and (multi-) critical points in a unified fashion, using the fact that all of them occur when 
the difference between phases with equal chemical potentials becomes infinitesimal. Because the truncatable models 
that we are considering are generally derived from mean- field approaches, we do not concern ourselves with critical 
point singularities, assuming instead that the free energy is a smooth function of all order parameters (the densities 
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in the system). Likewise, we do not discuss the subtle question of how, beyond mean- field theory, free energies can 
actually be defined in spinodal and unstable regions |^ . 

Let us first recap the general criteria for spinodals and critical points in multi-component systems. In order to treat 
the critera derived from the exact and moment free energies simultaneously, we use the common notation p for the 
vector of densities specifying the system: for the exact free energy, the components of p are the values p(cr); for the 
moment free energy, they are the reduced set pi. We write the corresponding vector of chemical potentials as 

m(p) = V/(p). 

This notation emphasizes that the chemical potentials are functions of the densities. The criterion for a spinodal at 
the parent phase p^^^ is then that there is an incipient instability direction 5p along which the chemical potentials do 
not change: 

{5p ■ V)p(p(o)) - {5p ■ V)V/(pW) = 0. (51) 

As the second form of the criterion shows, an equivalent statement is that the curvature of the free energy along the 
direction 5p vanishes. Eq. ( pT| ) can also be written as 

p(p(°)+6<5p)-p(pW) = 0(e2) (52) 

which is closely related to the critical point criterion that we describe next. 

Near a critical point, the parent p^°) coexists with another phase that is only slightly different; if, as we assume 
here, the free energy function is smooth, these two phases are separated - in p-space - by a "hypothetical phase" 
which has the same chemical potentials but is (locally) thermodynamically unstable. (This is geometrically obvious 
even in high dimensions; between any two minima of /(p) — p ■ p, at given p, there must lie a maximum or a saddle 
point, which is the required unstable "phase".) Now imagine connecting these three phases by a smooth curve in 
density space p(e). At the critical point, all three phases collapse, and the variation of the chemical potential around 
p(e = 0) = p^°^ must therefore obey 

p(p(e))-p(p(°)) ^0(e3) 

Similarly, if n phases coexist, we can connect them and the n — l unstable phases in between them by a curve p(e) [ p2[ . 
If we define an n-critical point as one where all these phases become simultaneously critical, we obtain the criterion 

p(p(6))-p(p("^) = 0(e2"-i). (53) 

This formulation was proposed by Brannock |^^; the cases n = 2 and n = 3 corres pon d to ordinary critical and 
tricritical points, respectively. The spinodal criterion (^2|) is of the same form as (|5^ ) if one chooses the curve 
p(e) = p(o) + e6p. 

In summary, we have that the phase p^^^ is a spinodal or n-critical point if there is a curve p(e) with p(e = 0) = p*-"-* 
such that 

Ap^p(p(e))-p(p(°)) = 0(e'). (54) 

where I — 2 for a spinodal, and / = 2n — 1 for an n-critical point. Brannock ]43| has shown that these criteria are 
equivalent to the determinant criteria introduced by Gibbs |^^. The above forms are more useful for us because 
they avoid having to define infinite-dimensional determinants (as otherwise required to handle the exact free energy 
/(p) of a polydisperse system, whose argument p denotes an infinite number of components, even in the truncatable 
case). They also show the analogy with the standard criteria for single-species systems (where p has only a single 
component) more clearly. 

We can now apply (g^ to the exact free energy ( |39| ) and show that the resulting criteria are identical to those 
obtained from the moment free energy. First, note that the curve p(e) = p{a; e) can always be chosen to lie within the 



family (40). This follows from the derivation of the criterion (p4|): The curve p(e) is defined as passing through I phases 
with equal chemical potentials, one of them being the parent p^^K As shown in (^9|), all these phases are therefore 
within the family (^). But then (^0|) implies that the condition ( ^4| ) that Afj,{a) must be zero to O(e') is equivalent 
to the same requirement for the differences Afii in the moment chemical potential. This proves that the conditions 
for spinodals and (multi-) critical points derived from the exact and moment free energies are equivalent. Intuitively, 
one can understand this as follows: Having shown that the spinodal/critical point conditions can be formulated solely 
in terms of density distributions p{a) within the family (^0|), it is sufficient to know the free energy of those density 
distributions. This is exactly the moment free energy. 
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Finally, we note that the exactness of these stability and critical point conditions holds not just for the parent 
p^'^\a), but for all phases within the family (|40|). This is clear, because any such phase could itself be chosen as 
the parent without changing the family; the moment free energy would change only by irrelevant terms linear in the 
moment densities. In general, one will not necessarily be interested in the properties of such "substitute" parents. As 
an exception, substitute parents that differ from the parent only by a change of the overall density are of interest: 
they lie on the dilution line of distributions from which cloud point and shadow curves are calculated. The dilution 
line is included in the family (40) if the overall density po (with corresponding weight function wo{a) = 1) is retained 
in the moment free energy; in the combinatorial derivation, this is automatically the case. 



B. Spinodals 



For completeness, we now give the explicit form ||11|,|12|,^ of the spinodal criterion (^1|) for truncatable systems; 
see also ]45| ] for an equivalent derivation using the combinatorial approach. Using (|4^) and abbreviating the matrix 
of second derivatives of the excess free energy as F, the spinodal condition dSl] ) becomes 



F + TM' 



0, 



{F + TM-'^)Sp = 0. 



(55) 



As before, the (nonzero) vector dp with components dpi gives the direction of the spinodal instability; the moment 
densities pi and pij that appear in F and M are to be evaluated for the parent distribution p*^"^ (a) being studied. 
Note that ( p5[ ) is valid for any p^'^'>{a); no specific assumptions about the parent were made in the derivation. We 
can thus simply drop the "(0)" superscript: for any phase with density distribution p{a), the point where ( [s^ ) first 
becomes zero, as external control parameters are varied, locates a spinodal instability. 

More convenient forms of (|5^) that avoid matrix inversions are obtained after multiplication by the second order 
moment matrix M (which is positive definite for linearly independent weight functions Wi{a) and therefore has 
nonzero determinant): 



Y = 1+ I3MF =0, (1 + f3MF)Sp = 0. 



(56) 



Here (3 = l/T in the standard notation. From our general statements in Sec. Ill A , the spinodal criterion derived 
from the exact free energy (^9|) must be identical to this; this is shown explicitly in App. Note that the spinodal 
condition depends only on the (first-order) moment densities pi and the second-order moment densities pij of the 
distribution p{a) (given by ( ^ ) and (]4^)); it is independent of any other of its properties. This simplification, which 
has been pointed out by a number of authors [Tl|JT^ , is particularly useful for the case of power-law moments (defined 
by weight functions Wi{a) = ct*): If the excess free energy only depends on the moments of order 0, 1. .. if — 1 of the 
density distribution, the spinodal condition involves only 2K — 1 moments [up to order 2{K — 1)]. 

The general discussion in Sec. [II A as well as the explicit calculation in App. ^ show that the spinodal instability 
direction lies within (or more precisely, is tangential to) the family (^) of density distributions. This fact has a simple 
geometrical interpretation: By construction, the family lies along a "valley" of the free energy surface (compare Fig. |lj). 
Away from the valley floor, any change in p{a) increases the free energy, corresponding to a positive curvature; the 
spinodal direction, for which the curvature vanishes, must therefore be along the valley floor. 



C. Critical points 

Next, we show the explicit form of the critical point criterion for truncatable systems. The general condition ( [s^ ) 
for an ordinary critical point {n = 2) was shown by Brannock p3| to be equivalent to 

(5p-V)V/(p("))=0, (,5p.V)V(p^°))=0 (57) 

The first part of this is simply the spinodal criterion, as expected. To evaluate the second part for the moment free 
energy (^8|), we need the third derivative of Sm w.r.t. the moment densities pi. Writing ( ^2| ) as 

^ dpjdpk 

and differentiating w.r.t. one of the moment densities, one finds after a little algebra 
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dpidpjdpk 

Iran 

with third-order moment densities pimn defined in the obvious way. Thus, evaluating (^7|) for the moment free 
energy (|8|), we find that critical points have to obey 

-T^v.iVjVk Pijk + SpiSpjSpk — = 0, Vi = ^(M"^),y Spj (58) 

ijk ap^apjOpk ^ 

in addition to the spinodal condition (|56|). As before, the criterion has to be evaluated for the parent phase p^'^\a) 
under consideration; but because p^^\(7) can be ch osen a rbitrarily, it applies to all density distributions p{<j). 



As expected from the general discussion in Sec. Ill A, the criterion (p8[) can also be derived from the exact free 
energy; an alternative form involving the spinodal determinant Y is given in App. Eq. ( |58| ) shows that the location 
of critical points depend only on the moment densities pi, pij and pijk [^l]j4^]. For a system with an excess free 
energy depending only on power-law moments up to order K — 1, the critical point condition thus involves power-law 
moments of the parent only up to order 3 (if — 1). 



D. Onset of phase coexistence: Cloud point and shadow 

So far in this section, we have shown that the moment free energy gives exact results for spinodals and (multi-) 
critical points. Now we consider the onset of phase coexistence, where (on varying the temperature, for example) a 
parent phase with density distribution ^("^((7) first starts to phase separate. As explained in Sec. |, this temperature 
together with the overall density p'"-* = p'j^^ of the parent defines a "cloud point" ; the density of the incipient daughter 
phase gives the "shadow". If the parent begins to coexist with p phases simultaneously (p = 2 at a triple point, for 
example), there will be p such shadows. 

At the onset of phase coexistence, one of the coexisting phases is by definition the parent p''^\(j); the lever rule does 
not yet play any role because the daughters p'-")(cr) {a — 1 . . .p) occupy an infinitesimal fraction of the total volume. 
It then follows from (^^ that all daughters lie within the family (|4^) of density distributions. As shown in (|5^), the 
condition for equality of chemical potentials /i(cr) between any two of the coexisting phases (parent and daughters) 
then becomes 

A/i(cr) = AfiiWtja) 

i 

and is satisfied if and only if the moment chemical potentials p,i are equal in all phases. Likewise, the exact and 
moment free e nerg ies give the same condition for equality of pressure 11 in all phases, because they yield identical 
expressions (^5[^8]) for 11. In summary, we see that the conditions for the onset of phase coexistence are identical for 
the exact and moment free energies; the moment free energy therefore gives exact cloud points and shadows. 



E. Phase coexistence beyond onset 



As stated in Sec. [I A , the moment free energy does not give exact results beyond the onset of phase coexistence, 
i.e., in the r egime where the coexisting phases occupy comparable fractions of the total system volume. As shown 
in Sec. [II A, the calculated phases will still be in exact thermal equilibrium; but the lever rule will now be violated 
for the "transverse" degrees of freedom of the density distributions. This is clear from (p^): In general, no linear 
combination of distributions from this family can match the parent p^^^ (cr) exactly. 

A more detailed understanding of the failure of the moment free energy beyond phase coexistence can be gained 
by comparing with the formal solution of the exact phase coexistence problem. Assume that the parent (a) has 
separated into p phases numbered by a = 1 . . .p. The condition (^), which follows from equality of the chemical 
potentials p{(j) in all phases, implies that we can write their density distributions p^'^\a) as 



P 



(a) = i?(a)exp(^^Af)u;.(a)^ 



(59) 



for some function R{a). If phase a occupies a fraction w^") of the system volume, particle conservation 
^^„(")p(")(cr) ^ p(o)(cr) then gives 
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exp (e, \^r^w,{a)) Y., exp (e, xf^w,{a)) 



If there are iiT moment densities pi, i — 1 . . . K, then this exact solution is parameterized by {p— l)(-fir + 1) independent 
parameters: {p— 1)K parameters A^"'' (noting that the Aj of one phase can be fixed arbitrarily), and p—l parameters 
yia) (^noting that one phase volume is fixed by the constraint J2a = Comparing ( |59| ) with (p9[), one sees that 
the K quantities 

TAf)+Mf^ {^ = 1...K) 

must be the same in all phases a; the same is true for the pressures n("\ and this gives the required total number of 
(p — 1){K + 1) constraints. This formally defines the exact solution of the phase coexistence problem for truncatable 
systems. Its practical value is limited by the difhcultics of finding the solution numerically, as pointed out in Sec. |; 
see also Sec. |V|. There is also no interpretation of the result in terms of a free energy depending on a small number 
of densities. 

Nevertheless, eq. ( |60| ) is useful for a comparison with the solution provided by the moment free energy method. 
From (59), the exact coexisting phases are all members of a family of density distributions of the general form (|^); 



this is simply a consequence of the requirement of equal chemical potentials ([4 91) . This "exact coexistence family" has 



an "effective" prior |47| i?(cr) ^ p'^'^\a) and is therefore different from the "original" family (pO|). The moment free 
energy only gives us access to distributions from the original family, not the exact coexistence family, and therefore 
cannot yield exact solutions for phase coexistence (beyond its onset). Fig. |^ illustrates this point explicitly in the 
simplified context of a bidisperse system. 



It is now easy to see, however, that - as stated in Sec. [I A - the exact solution can be approached to arbitrary 



precision by including extra moment densities in the moment free energy. (This leaves the exactness of spinodals, 
critical points, cloud-points and shadows unaffected, because none of our arguments excluded a null dependence of 
/ on certain of the pi.) Indeed, by adding further moment densities one can indefinitely extend the family (^) of 
density distributions, thereby approaching with increasing precision the actual distributions in all phases present; this 
yields phase diagrams of ever-refined accuracy. 

The roles of the "original" moment densities (those appearing in the excess free energy) and the extra ones are 
quite different, however. To see this, note from (E3) that equality among phases of the moment chemical potentials 
implies that, for the extra moments only, the corresponding Lagrange multipliers must themselves be equal in all 
coexisting phases. (This is because there is, by construction, no excess part to the "extra" chemical potentials.) We 
can therefore drop the phase index on these Lagrange multipliers and write the density distribution in phase a as 



p(")(a) = p(°)(a)exp V X,w,{a) exp V x't^w,{a) (61) 





Comparing with (0) and (^9|), we see that the extra Lagrange multipliers can be thought of as providing a "flexible 
prior" R{a), allowing a better approximation to the effective prior R{cr) required by particle conservation. The 
"tuning" of the prior with these extra Lagrange multipliers also has an important effect on the number of coexisting 



phases that can be found by the moment method (see Sec. IIIG) 



F. Global and local stability 



Most numerical algorithms for phase coexistence calculations, including ours, initially proceed by finding a solution 
to the equilibrium conditions of equal chemical potentials and pressures in all phases, rather than by a direct mini- 
mization of the total free energy of the system. It is then crucial to verify whether this solution is stable, both locally 
{i.e., with respect to small fluctuations in the compositions of the phases) and globally with respect to splitting into 
a larger number (or different) phases (see Fig. ^). Global stability is a particularly important issue in our context 
because there is in principle no limit on the number of coexisting phases in a polydisperse system. 

A useful tool for stability calculations is the "tangent plane distance" [Q. Let us first define this generically 
for a system with (a vector of) densi ties p, free energy density f{p) and chemical potentials fJ,{p) = V/(p); this 
notation is the same as in Sec. [II A. Assume we have found a candidate "phase split", that is a collection of p 
phases p^"-* (a — 1 . . .p) which satisfy the phase equilibrium conditions of equal chemical potentials p — p'^"^ and 
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FIG. 2. This figure continues the example of Fig. which is a Flory-Huggins polymer+solvent model with two chain lengths 
present. Now, however, we allow for a change of the interaction parameter x (equivalent to changing temperature), (a) The 
moment free energy at the cloud point, i.e., the value of x for which the given parent first begins to phase separate, (b) The 
tieline at the cloud point is drawn in the (n, r2) plane of all density pairs (thin solid line). It is exact and connects the parent 
(solid circle) with another member of the family (|4o|); the family is indicated by the thick solid line. The dotted line is the 
"transverse" line of constant pi (= p'°') passing through the parent, (c) The moment free energy for a higher value of x- (d) 
The corresponding tieline (thin dashed line) in the (ri, r2) plane. As before, this connects two members of the parent's family 
(now the thick dashed line), and the two phases so found are in stable thermodynamic equilibrium with each other. Because 
they now occupy comparable fractions of the system volume, however, the lever rule is violated: the parent does not lie on 
the tieline. But the lever rule violation occurs only along the transverse direction (dotted line): the fractional phase volumes 
calculated from the moment free energy result in a total density pair (indicated by the empty circle) which has the same moment 
density p as the parent. Finally, the exact tieline at this x is also shown (thin solid line). As required, it passes through the 
parent. But its endpoints now connect members of a different family (thick solid line), which derives from an "effective parent" 
(or effective prior) (fi, f2), which may be chosen anywhere along the thick solid line (compare (^^). Note that if a single extra 
moment density were added in this scenario, the resulting maximum entropy family would in fact cover all possible density 
pairs (ri, r2) and the corresponding two-moment free energy would therefore give exact results in all situations. For a truly 
polydisperse (rather than bidisperse) system, in which the density pairs of this example become density distributions p{(j), each 
added moment density allows the accuracy of the calculation to be increased. 
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FIG. 3. Illustration of local and global stability of solutions of the phase equilibrium conditions. For simplicity, we consider 
the case of a monodisperse system, characterized by its free energy density / (bold line) as a function of the density p. The 
phase split for a parent with density p'°' (indicated by the circle) is to be calculated. Three double tangents to /(p) are shown, 
predicting different two-phase splits of the parent. All obey (by the double tangent property) the conditions of equal chemical 
potential and pressure in the two phases. But (1) is locally unstable: a small fluctuation (indicated by the arrow) of the density 
of the more dilute coexisting phase lowers the total free energy. (2) is locally stable, but globally unstable: the phase split (3) 
gives a lower total free energy. (3) is globally (and therefore locally) stable. 

pressures n = n("). They occupy fractions t;'"-* of the total volume, and the overall density distribution is thus 
Ptot = p'"''- Note that we do not yet assume that the lever rule is satisfied, i.e., that p^ot = P^^^'i explained 

above, this equality will generally not hold for phase splits calculated from the moment free energy. We now define 
"global tangent plane (TP) stability" for such a phase split as the property that there is no other phase split (i.e., 
no other tangent plane) that gives a lower total free energy for the same overall density distribution p^^^. In more 
intuitive language, this means that if we were to put the phases p'"-* into contact with each other, the resulting system 
would be thermodynamically stable; neither the composition nor the number of phases would change over time. Note, 
however, that since in general p^^^ = ■(;(") is not equal to p'°', the phase split that is globally TP-stable need 
not accurately reflect the number and composition of phases into which the parent p^'^-' would actually split under 
the chosen thermodynamic conditions. 

To define the tangent plane distance (TPD), note first that by virtue of the coexistence conditions, all phases p^"^ 
lie on a tangent plane to the free energy surface. Points (p, /) on this tangent plane obey the equation / — p-p+II = 0, 
with p and 11 the chemical potentials and pressure common to all phases. For a generic phase with density distribution 
p and free energy /(p), the same expression will have a nonzero value which measures how much "below" or "above" 
the tangent plane it lies. This defines the TPD 

i(p) = /(p)-p(").p + n("). (62) 

Here we have added the superscript a to emphasize that the chemical potential and osmotic pressure used in the 
calculation of the TPD are those of the calculated phase equilibrium (and hence of any of the participating phases 
a), rather than those of the test phase p. It is then clear intuitively - and can be shown more formally pC| - that 
the calculated phase coexistence is globally TP-stable if the TPD is non-negative everywhere. Geometrically, this 
simply means that no part of the free energy surface must protrude beneath the tangent plane; otherwise the total 
free energy of the system could be lowered by constructing a new tangent plane that touches the protruding piece. 
Global TP-stability of course encompasses local stability; the latter simply corresponds to the requirement that the 
TPD be a local minimum (with the value t = 0) at each of the phases in the candidate solution. 

When verifying global TP-stability, it is obviously sufficient to check the value of the TPD at all of its stationary 
points. By differentiating (^) w.r.t. p, one sees that at these points, the chemical potentials p(p) are the same as in 
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the calculated coexisting phases. For our truncatable polydisperse systems, it then follows from eq. ( ^9| ) that we only 
need to consider the TPD for test phases p{<7) which are in the same family (0) as the calculated coexisting phases. 
This is a crucial point: Even though global TP-stability is a statement about stability in the infinite-dimensional 
space of density distributions p(cr), it can be checked by only considering density distributions from a ii'-dimensional 
family. In the (generally hypothetical) case of an exactly calculated phase split, this family would be ( pOj) , with the 
effective prior i?(iT) For a phase split calculated from the moment free energy, it is the family ( pio|) with the 

parent as prior, with any Lagrange multipliers for extra moments fixed to their values in the calculated coexisting 
phases. The TPD of such a test phase is (using (^9|) for the exact free energy and (|4^ ) for the chemical potentials 
pS°'\ct) of the coexisting phases) 

t[p{a)] = fW)] - I dap'^"\a)p{a) + n'") 



= /(p.) + j dap{a) |t [Inp(a) - 1] - T\np^^\a) ~ ^ ~^^t^w.{o)^ + H^") 



f{p.)+T jdapia) p^(o)(^) 

/m(pO-EMi"^'^^ + n'"'- 



A comparison with ( |62| ) shows that this is identical to the TPD that one would derive from the moment free energy 
alone. Here again, the underlying polydisperse nature of the problem can therefore be disregarded once the moment 
free energy has been obtained. In summary, one can determine whether a phase split calculated from the moment free 
energy is globally TP-stable (which means that phases of the predicted volumes and compositions, would, if placed 
in contact, indeed coexist) using only the TPD derived from the moment free energy; the same is trivially true of the 
weaker requirement of local stability |^ . 

Recall, however, that the overall density distribution ptot(o') for a phase split found from the moment free energy only 

has the same moment densities pf^ as the parent [a) , but differs in other details (the transverse degrees of freedom) . 
Global TP-stability thus guarantees that such an approximate phase split is thermodynamically stable, but does not 
imply that it is identical (in either number or composition of phases) to the exact one. Nevertheless, progressively 
increasing the number of extra moment densities in the moment free energy will make ptot (c) a progressively better 
approximation to p*^"^ (ct), and so the exact phase split, with the correct number of phases, must eventually be recovered 
to arbitrary accuracy; see Fig. ^ for an illustration. 



G. Geometry in density distribution space 



In the above discussion of the properties of the moment free energy, we have focused on obtaining the phase behavior 
of a system with a given parent distribution p^'^^ (cr) . This is the point of view most relevant for practical applications 
of the method, and the rest of this Section is not essential for understanding such applications. Nevertheless, from a 
theoretical angle, it is also interesting to consider the global geometry of the space of all density distributions p{<j), 
without reference to a specific parent: the exact free energy ( |39| ) of the underlying (truncatable) model induces two- 
phase tie lines and multi-phase coexistence regions in this space, and one is led to ask how the moment free energy 
encodes these properties. 

As pointed out in Sec. II A, the definition of the moment free energy depends on a prior R{a) and represents the 
properties of systems with density distributions p{f7) in the corresponding maximum entropy family (^. Instead 
of identifying R{a) — p'-*'-' (cr) , we now allow a general prior R{a). Conceptually, it then makes sense to associate 
the moment free energy with the family (0) rather than the specific prior. This is because any density distribution 
from ^ can be chosen as prior, without changing the moment free energy (apart from irrelevant linear terms in the 
moments p^), or the identity of the remaining family members. The construction of the moment free energy thus 
partitions the space of all p(cr) into (an infinite number of) families (^; different families give different moment free 
energies that describe the th ermod ynamics "within the family" . This procedure gives meaningful results because, 
as shown at the end of Sec. Ill A, coexisting phases are always members of the same family. By considering the 
entire ensemble of families (^) and their corresponding moment free energies, one can thus in principle recover the 
exact geometry of the density distribution space and the phase coexistences within it. Note that this includes regions 
with more than K + 1 phases, even when the excess (and thus moment) free energy depends on only K moment 
densities; see Fig. |[ Consistent with Gibbs' phase rule, however, the families for which this occurs are exceptional: 
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FIG. 4. Example of structure of the space of density distributions. We continue the example of Figs. ^ and |2| and consider 
a bidisperse Flory-Huggins system with a density pair (n, r^). As before, we take the excess free energy of the system to 
depend on a single moment density, so that the families (|7|) are one- dimensional {i.e., appear as curves in the graph). In the 
case sketched here, there is a three-phase region bordered by two-phase regions delineated by dotted lines. The bold (solid, 
dashed, grey) curves show three families (Q) of density pairs; the whole space is partitioned (or "foliated") into an infinite 
number of such families. All tielines must begin and end within the same family; the thin grey tielines are shown as examples. 
The moment free energy (without use of extra moments - or in this case, with a single moment density retained) can only 
access systems from the family passing through the parent (filled circle); it therefore predicts that the parent will separate 
into the two phases within this family that are connected by a tieline (thin dashed line). This tieline is thermodynamically 
"real": the two phases at its ends, if put into contact, would remain in stable coexistence. But it does not give the exact phase 
separation for the given parent because of the lever rule violations allowed by the moment approach (the parent does not lie 
on the tieline). Retaining one extra moment density in the moment free energy gives, in this simple bidisperse case, the exact 
result: the parent separates into the three phases at the corners of the three-phase triangle, which lie on a family (bold solid 
curve) that does not contain the parent. 
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they occupy submanifolds (of measure zero) in the space of all families. Accordingly, to find the families on which 
such "super-Gibbs" multi-phase coexistences occur, the corresponding prior must be very carefully tuned. Indeed, 
the probability of finding such priors "accidentally" , without either adding extra moments to the moment free energy 
description or solving the exact phase coexistence problem, is zero. Generically, one needs to retain at least n moment 
densities in the moment free energy to find n + 1 phases in coexistence. From (|6T|), this corresponds to having n — K 
parameters available for tuning the location of the family (or equivalently, its prior) in density distribution space. 
But note that, once the extra moments have been introduced, the required tuning need not be done "by hand": it is 
achieved implicitly by requiring the lever rule to be obeyed not only for the K original moment densities, but also for 
the n — K extra ones. (The solution of the phase equilibrium conditions thus effectively proceeds in an n-diniensional 
space; but once a solution has been found, its (global T P-)sta bility can still be checked by computing the TPD only 



within a if-dimensional family of distributions. See Sec. [II F.) 



As pointed out above, phase splits calculated from the moment free energy allow violations of the lever rule. In the 
global view, this fact also has a simple geometric interpretation: the families (^ are (generically) curved. In other 
words, for two distributions /9'^^-'((t) and /9^^-'(cr) from the same family, the straight line p{a) — ep^^''{a) + (1 — e)p^^\a) 
connecting them lies outside the family. More generally, if p coexisting phases p^°'\a) have been identified from the 
moment free energy, then the overall density distribution /Otot(o') = v'^"^ p^"^ {a) is different from the parent /^'"■'(cr); 
geometrically, it lies on the hyperplanc that passes through the phases p^°'\a). 

IV. PRACTICAL IMPLEMENTATION OF THE MOMENT METHOD 

The application of the moment fr ee e nergy method to the calculation of spinodal and critical points is straighforward 
using conditions (|5^), and ( |58| ) or (|D3|), respectively, and is further illustrated in Sec. ^ below. We therefore focus in 
this Section on phase coexistence calculations. 

Recall that in the moment approach, each phase a is parameterized by Lagrange multipliers A^"'' for the original 
moments (the ones appearing in the excess free energy of the system) and the fraction u^"^ of system volume that 
it occupies. If extra moments are used, there is one additional Lagrange multiplier for each of them; these are 
common to all phases. These parameters have to be chosen such that the pressure (^ ) and the moment chemical 

are equal in all phases. Furthermore, the (fractional) phase volumes w^"^ have to sum to 
las to be satisfied for all moments (both original and extra): 

\ ^ (a) (a) (0) 



potentials pi given by 
one, and the lever rule 



In a system with K original moments that is being studied using an n-moment free energy {i.e., with n — K extra 
moments), and for p coexisting phases, one has p{K + l) + n — K = {p— 1){K + 1) + n + 1 parameters and as many 
equations. Starting from a suitable initial guess, these can, in principle, be solved by a standard algorithm such as 
Newton- Raphson EG]. Generating an initial point from which such an algorithm will converge, however, is a nontrivial 
problem, especially when more than two phases coexist. 

To simplify this task, we work with a continuous control parameter such as temperature, density po of the parent 
phase, or interaction parameter x for polymers. Taking the latter case as an example, we start the calculation at a 
small value of x where we are sure to be in a single-phase region; thus the parent phase under consideration is stable. 
The basic strategy is then to increment x ^'iid detect potential new phases as we go along. At each step, the phase 
equilibrium conditions are solved by Newton-Raphson for the current number of phases. Then we check for local 
stability of the solution by calculating the Hessian of the TPD around each of the coexisting phases and verifying 
that it is positive definite. If an instability {i.e., a negative eigenvalue of the Hessian) is found, we search for local 
minima starting from points displaced either way along the instability direction. If two new local minima are found 
in this way, we add them to list of phases and delete the unstable phase; if only one new local minimum is uncovered. 



we add this but retain the old phase. Finally, we check for global (TP-)stability. As explained in Sec. IIIF, this 
involves scanning all "test" phases from the same maximum entropy family for possible negative values of the TPD. 
However, the test phase will have the same extra Lagrange multipliers, and is therefore parameterized in terms of 
the Xi for the original moments only. The extra Lagrange multipliers are held fixed, so their number is irrelevant for 
the computational cost of the stability check. Put differently, within our algorithm the precision of the coexistence 
curves should depend on the total number of moments retained, n, rather than the number of moments in the excess 
free energy, K. Computational effort, on the other hand, is dominated by the global stability check, and thus mainly 
sensitive to K (the dimension of the space to be searched for new phases), rather than n. 

This is a substantial efficiency gain, but despite it, an exhaustive search for local TPD minima over the K- 
dimensional space of original Lagrange multipliers is unrealistic except for AT = 1. We use instead a Monte-Carlo- type 
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algorithm to sample the TPD at representative points. If any negative values of the TPD are encountered, we 
choose the smallest such value, find the nearest local minimum of the TPD and add this point to the list of phases. 
If any new phases have been found during the (local and global) stability checks, we assign each of them a default 
phase volume [v = 0.1, say), reduce the phase volumes of the old phases accordingly, update the number of phases 
and loop back to the Newton-Raphson solution of the phase equilibrium conditions, using the current list of phases 
as an initial guess. This process is repeated until the calculated solution is found to be stable. Our implementation 
of this basic scheme also contains some additional elements (such as an adaptive choice of the stepsize for the control 
parameter x, and checks for very small phase volumes) which are useful near points where new phases appear or old 
ones vanish. 

Comparing our own approach, which is based on the moment free energy with extra moment densities as just 
outlined, to (|60|), which is the exact solution for truncatable systems, we see that the former actually uses more 
parameters (essentially one phase-independent Lagrange multiplier per extra moment) to represent the solution. At 
first sight, this may appear counter-productive. However, the Lagrange multipliers A^"^ and the (fractional) phase 
volumes w'"-* are much less strongly coupled in the moment free energy approach. The phase volumes are only 
determined by the lever rule, while in the exact solution they "feed back" into the effective prior R{a) and therefore 
into the equilibrium conditions of equal chemical potentials and pressures. In the examples studied below, we have 
found that the numerical advantages of this decoupling (in terms of stability, robustness and convergence of our 
algorithm) can easily outweigh the larger number of parameters that it requires. 

As with any numerical algorithm, it is important to develop robust criteria by which the convergence of the solution 
can be judged. This has a novel aspect, in the moment method, since violations of the lever rule are allowed: alongside 
normal numerical convergence criteria one needs a method for decidin g wh ether the effect of these on the predicted 
phase behaviour is significant. We develop appropriate criteria in Sec. VB in the context of a specific example. The 
basic ideas is that since any state of phase coexistence predicted by the method represents the exact behaviour of 
some parent, then so long as this parent is close enough to the true one, the predicted behaviour will lie within the 
range of uncertainty that arises anyway, from not knowing the true parent to arbitrary experimental precision. 



V. EXAMPLES 



We now illustrate how the moment method is applied and demonstrate its usefulness for several examples. The 
first two (Flory-Huggins theory for Icngth-polydisperse homopolymers and dense chemically polydisperse copolymers, 
respectively) contain only a single moment density in the excess free energy and are therefore particularly simple to 
analyse and visualize. In the third example (chemically polydisperse copolymers in a polymeric solvent), the excess 
free energy depends on two moment densities and this will gives us the opportunity to discuss the appearance of more 
complex phenomena such as tricritical points. 



A. Homopolymers with length polydispersity 

Let us start with the simple but well studied example of polydisperse Flory-Huggins theory ||5^]. One considers a 
system of homopolymers with a distribution of chain lengths; the polydisperse feature a is simply the chain length 
L, i.e., the number of monomers in each chain. (We treat this as a continuous variable.) The density distribution 
p{L) then gives the number density of chains as a function of L. We choose the segment volume as our unit of 
volume, making p(L) dimensionless. The volume fraction occupied by the polymer is then simply the first moment 
of this distribution, 4> = pi = JdL L p{L). Within Flory-Huggins theory, the free energy density is (in units such that 
fceT = 1) 

/ = j dLp{L)[\np{L)-l] + {l-p^)\n{l^p{) + xPi{l-pi) (63) 

where the Flory ^-parameter plays essentially the role of an inverse temperature. Before analysing this further, we 
note that for a bidisperse system with chain lengths Li and L2 and number densities ri and r2, the corresponding 
expression would be 

/ = ri(lnri - 1) + r2(lnr2 - 1) + (1 - pi) ln(l - pi) + xpi(l - Pi), Pi = ^i^i + L2r2 (64) 

This free energy, with Li — 10 and L2 = 20, was used to generate the examples shown in Figs. |^ and ^. 
Returning now to (|63|), we note first that the excess free energy 
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/ = (1 - Pi) ln(l - pi) + xpi(l - pi) 



(65) 



depends only on the moment density pi. If no extra moments are used, the moment free energy is therefore a function 
of a single density variable, pi. Let us work out its construction explicitly for the case of a parent phase with a 
Schulz-distribution of lengths, given by 



1 



r(a)6'' 



(66) 



Here a is a parameter that determines how broad or peaked the distribution is; it is conventionally denoted by a, 
but we choose a different notation here to prevent confusion with the phase index used in the general discussion so 
far. The chain number density of the parent is p^^\ and the normalized first moment m^*^^ — pf''^ / p'j^'^ is simply the 
(number) average chain length 



(0) 

r - (0) _ Pi . 
Ln ^ m\ ' = = ab. 



(67) 



The parent is thus parameterized in terms of a, pl^^ and p''^'' (or b). The density distributions in the family (40) are 
given by 
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p(i) =p(")(L)e/^^ =p(")^L-ie-(^-^-^^)^ 



(68) 



and the zeroth and first moment densities are easily worked out, for members of this family, to be 

a+l 



(0) 
po = Po 
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(0) 

Pi = Pi 
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1 - b\l 



Because the family only has a single parameter Ai, these two moment densities are of course related to one another: 

-a/(a+l) 
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(69) 



Now we turn to the moment entropy, given by (|3q): Sm — po — Aipi, and to be considered as a function of pi. After 



a little algebra (and eliminating b in favor of Pi) one finds 
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(70) 



The last term is linear in pi and can be disregarded for the calculation of phase equilibria. This then gives the 
following simple result for the moment free energy 



/„, = -(a + l)cp?/("+^^+/. 



(71) 



(In conventional polymer notation, this result would read /,n = —{a + l)c0"/("+^) + /.) From this we can now obtain 
the spinodal condition, for example, which identifies the value of x where the parent becomes unstable. The general 
criterion ( ^l|) simplifies in our case of a single moment density to 
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(72) 



The same condition for a phase with general density distribution p{L) - rather than our specific Schulz parent p^^^ (L) 
- follows from (^) as 



1 „ 1 
2x + — 

1 - Pi P2 



(73) 



As expected, this becomes equivalent to (^ for parents of the Schulz form (|66[), which obey p'^^ — p'^\{a + 1)5^ 



and eq. (|67|). Note that in standard polymer notation, eq. ([7 3D would be written as 
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(a) X=0.55 



(b) 5(=0.585 



(c) X=0.62 




0.05 0.1 0.15 0.2 0.05 0.1 0.15 0.2 0.05 0.1 0.15 0.2 

pi pi pi 

FIG. 5. Examples of moment free energy ( j7l| ) for Flory-Huggins theory of length-polydisperse polymers, with one moment 
density, pi, retained. The parent is of the Schulz form (|66|), with pj"' = 0.03, I/jv = 100 (hence Pg"' = Pi'' I^n = 3 • 10"*) and 
a = 2 (hence Lw = 150); the point pi = p^^^ is marked by the filled circles. In (a), the value of x = 0.55 is sufficiently small 
for the parent to be stable: The moment free energy is convex, (b) shows the cloud point, x ~ 0.585, where the parent lies on 
one endpoint of a double tangent; the other endpoint gives the polymer volume fraction pi in the shadow phase. Increasing x 
further, the parent eventually becomes spinodally unstable (x ~ 0.62, plot (c)). Note that for better visualization, linear terms 
have been added to all free energies to make the tangent at the parent coincide with the horizontal axis. 



1 — '''' Ly/(f) ^ ^ 

where — P2I P\ is the weight average chain length. 

The critical point condition is obtained similarly. For a single density variable, eq. (m) reduces to 



dpi 



aia + 2) ^ ^^^^ 



pi=p\ 

For a general distribution, the critical point criterion (Esh becomes instead 



pI dpi 



= 



and inserting the Flory-Huggins excess free energy (|6^) gives 

1 _ P3 
{l^Pi? P\ 



(76) 



Again, this simplifies to (|75| ) for Schulz parents, for which pg^"* = p^^^ a{a + l)(a + 2)6'^. The traditional statement 
of d?^), in terms of </> = pi, Ly/ and the Z-average chain length Lz = P3/P2 reads [ p2[ 

(WF"lf0^ = °- 

Fig. ^ shows the moment free energy ( |7l| ) for a given parent and different values of x- The value of x for which 

there is a double tangent at the point pi = p^"^ representing the parent gives the cloud point; the second tangency 
point gives the polymer volume fraction pi in the coexisting phase, hence the shadow point. Repeating this procedure 
for different values of p^°^ (while maintaining Lpf = Pi°VPo°'' ^^'^ constant), one obtains the full CPC and shadow 
curve for a given normalized parent length distribution n(L) — p^'^\L)/pl^\ 

We now consider the properties of the moment free energy with the (chain number) density po retained as an 
extra moment. This provides additional geometrical insight into the properties of polydisperse chains (while for the 
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numerical determination of the CPC and shadow curve, the above one-moment free energy is preferable). To construct 
the two-moment free energy, we proceed as before. The family (^o|) is now 



p(L)=p(°)(i)e 



with zeroth and first moment densities 
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The relations can be inverted to express the Lagrange multipliers Aq and Ai in terms of the moment densities as 

Ao = (a + 1) In - a In 

Po Pi 



Ai = 

and one obtains the moment entropy 

Sm = Po - Aopo - Aipi = -(a 
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(79) 



As expected, this reduces to ([70|) for systems with density distributions from our earlier one-parameter family (|6q), 
where po can be expressed as a function of pi according to (69). With both po a-nd pi retained as moment densities 
in the moment free energy, a number of linear terms in (|79|) can be dropped, giving the final result 



./m = (a + l)po Inpo - apo hipi + /. 



(80) 



Note that the dependence on the parent distribution is now only through a. This can be understood from the general 
discussion in Sec. Ill Gj : The family ( [T^ ) of density distributions now contains all Schulz-distributions with the given 
a, and the moment free energy is insensitive to which member of this family (specified by p^^^ and pf^) is used as the 
parent. 

The result (|80| ) can of course also be obtained via the combinatorial method, as follows. The normalized Schulz 
parent distribution is given by n(°)(L) = [r(a)6"]~^L°~^e~''^. The cumulant generating function is thus h{\i) = 
—a ln(l — fcAi), giving mi = dh/dXi — ab/ (1 — bXi). Dropping constants and linear terms in the average chain length 
TOi (which do not affect the phase behavior), the generalized entropy of mixing per particle, Smix = h — Aimi then 
becomes Smix = alnmi. Adding the ideal gas term and the excess part of the free energy, we thus find the moment 
free energy density 



/,„ = Po Inpo - apo ln(pi/po) + / 



(81) 



in agreement with 



it is a simple two component free energy. In the conventional notation, eq. (81) would read 
/„! = plnp - ap\n{(f)/p) + {1 - (p) ln(l - (j)) + x0(l - 0). 



The spinodal curve (SC) and critical point (CP) condition may now be calculated from eq. ( pl|) , using either the 
methods outlined in Sec. pll A or the more traditional determinant conditions |Q (see also App. y). One finds 
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0, (SC) 
(CP) 



(82) 
(83) 



Evaluated at the parent, these conditions are identical to eqs. ( |7^ ) and ( k5|) - which were derived from the one- 
moment free energy - as they must be. But, as explained at the end of Sec. Ill A, they also hold more generally for 
all systems with Schulz density distributions with the given value of a; see the remarks after eq. (80). Of particular 
interest among these systems are those that differ from the parent only in their overall density while having the same 
normalized length distribution n{L) — n^^\L). They form the "dilution line"; their number averaged chain length 



26 




FIG. 6. Conventional two component phase behavior in polydisperse Flory-Huggins theory, shown in the {pi,po) plane for 
three values of x- As in Fig. ^ the parent has Ln ~ 100 and Lw = 150 (hence a — 2). Along the j/-axis, we plot Ljvpo 
rather than po so that the dilution line pi = Ljvpo, shown as the thick solid line in (a-c), is simply along the diagonal. With 
X considered as an additional variable, the dilution line constraint defines a plane {pi — LnPo,x)- The last plot, (d), shows 
the cut by this plane through the phase behavior in (a-c); the solid line is the cloud point curve, and the dashed line is the 
spinodal stability condition. 



satisfies uii = m^"-' = J dLLn^^\L), which translates to pi = L^po (where Ljsi = mf^ = ah is, as before, the number 
average chain length of the parent). Inserting this "dilution line constraint" into the above equations, we obtain as 
the spinodal and critical point conditions for systems on the dilution line 

^ 2x+^-^=0, (SC) (84) 



Pi g+ILnPi 
1 a{a + 2) 1 



xo / 9=0. (CP) (85) 

(l-pi)2 (a + 1)2 Lwp2 V 7 \ ) 

Again, agreement with the general results (^4|) and (^^ is easily verified by noting that for Schulz length distributions, 
Lw — P2/P1 = (a + l)b = L]^{a + l)/a and Lz = Pzj P2 = (a + 2)6 = Ljq{a + 2)/a for the weight and Z-average chain 
lengths, respectively. 

The above discussion focussed on spinodals, critical points, and cloud and shadow curves, all of which are found 
exactly by the moment approach. (As usual, "exact" is subject to the assumed validity of the original, truncatable 
free energy expression (|65[).) Within the moment approach, we have also calculated the full phase behavior in the 
(piiPo) plane. This was done numerically from the two-moment free energy (^l]), facilitated by a partial Legendre 
transformation from po to the associated chemical potentialuo, which can be performed analytically in this case (and 
effectively brings one back to the one-moment free energy ([7l|)). Fig. ^(a)-(c) shows binodal curves, tielines, spinodal 
curves and critical points in the (pi,po) plane, for three values of x- These plots give an appealing geometric insight 
into the problem of phase separation in this system. For instance, the slope of the tielines indicates a size partitioning 
effect: the more dense phases are enriched in long chains. 

The heavy line in Fig. |^(a)-(c) is the dilution line constraint rrii = wl'^^ or p\ — L^pQ. As discussed already, we are 
interested mainly in systems whose mean composition lies on this line. (With x as an extra variable, this line becomes 
a plane (pi = LnPo, x) in the space {pi, po, x)-) Not all of the phase behavior shown in Fig. ^ is then accessible. The 
extremities of phase separation on the pi — Ljvpo line are points where phase separation just starts to occur, and the 
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FIG. 7. The phase separation region, and compositions of the shadow phases corresponding to the two cloud points, for 
X = 0.62 (compare Fig. n(c)). 




FIG. 8. Phase behavior in (pi — LnPo,x) plane showing (a) pi-projected shadow curve, and (b) po-projected shadow curve. 



locus of these points in the (pi = LatPo: x) plane defines the cloud curve, which bounds the phase separation region. 
This is shown in Fig. ||(d). This plot also includes the spinodal stability curve and the critical point from eqs. (jsj) 
and (p5[). The spinodal curve and the cloud curve touch at the critical point, which no longer lies at the minimum 
of either. This distorted behavior is a well known feature of polydisperse systems. Here it is seen to be due to the 
way that regular phase behavior in (pi,/9o,x) space is cut through by the dilution line constraint. Note that only 
two moments are needed to understand this qualitative effect, and that Fig. ^ gives direct geometrical insight into it, 
which would be very hard to extract from any exact solution based on (j5^) and (|60[). 

The compositions of the phases that just start to appear as the phase separation region is entered do not in general 
lie on the dilution line pi = L^pq. These phases are the "shadows"; they lie at the other ends of the cloud point 
ticlines in Fig. ^(a)-(c). This is illustrated for x — 0-62 in Fig. 0. The shadow phase compositions may be projected 
onto the (pi — Ljvpo, x) plane to give shadow curves, but the projection is not unique. The shadow curve obtained 
by ignoring the value of po and projecting onto the pi-axis, for example, is shown in Fig. |^(a). That obtained by 
ignoring the value of pi and projecting onto the po-axis is shown in Fig. |8|(b). The different shapes are due to the fact 
that the locus of shadow phase compositions in general only intersects the (pi — L^po^x) plane at a single point — the 
critical point. 

In conclusion, we emphasize again that all the curves shown in Fig. |^(d) and Fig. ^ are exact, even though they 
have been constructed from the phase behavior shown in Fig. ^(a)-(c) which is itself a projection of the true phase 
behavior of the fully polydisperse system. 



B. Copolymer with chemical polydispersity 

Next, we apply the moment free energy formalism to random AB-copolymers in the melt or in solution For 
simplicity, we assume that the copolymer chain lengths are monodisperse, so that each chain contains an identical 
number L of monomers. Even with this simplification, the problem is quite challenging because, as we shall see, a 
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given parent can split into an arbitrarily large number of phases as the interaction strength is increased. 

For this problem we choose as our polydisperse feature, cr, the difference between the fractions of A and B-monomers 
on a chain. This takes values in the range cr = —1 ... 1, with the extreme values ±1 corresponding to pure A and 
B chains, respectively. We now measure volumes in units of the chain volume (which is L times larger than the 
monomeric volume used in the previous example). The total density po — Jd<T p{a) is then just the copolymer volume 
fraction (p. (Note that for the case of length polydispersity, (j) was related to the first moment of the length distribution; 
here it is the zeroth moment of the distribution of chain compositions cr.) Within Flory-Huggins theory, the free energy 
density (in units of k^T per chain volume) is 

/ = jdap{a)[\np{a)-l] + -^{l-po)Hl-po) + Lxp^{l-po)-Lj('pl (86) 

The interaction between A and B monomers enters through pi = Jdacrp{cr). In (|8^), we have specialized to the case 
where the interactions of the A and B monomers with each other and with the solvent are symmetric in A and B; 
otherwise, there would be an additional term proportional to poPi- We have also allowed the solvent to be a polymer 
of length Lg; the case of a monomeric solvent is recovered for Lg — 1. In the following, we absorb factors of L into 
the interaction parameters by setting 

X = Lx, x' = Lx' 
Defining also r — LjL^^ the excess free energy takes the form 

/ - r(l - po) ln(l - po) + XPo(l - Po) - x'p?- (87) 
The moment free energy (p8[), with only po and p\ retained, is then 

/m Aopo + Aipi - Po + '"(1 - Po) ln(l - Po) + XPo(l - Po) - x'Pi- (88) 
If extra moment densities with weight functions Wi{a) are included, these simply add a term XiPi each. 



1. Copolymer without solvent 

We first consider the special case of a copolymer melt, for which the overall density (i.e., copolymer volume fraction) 
is constrained to take everywhere its maximum value po = 1. The free energy (^) then simplifies to 

/m = Ao + Aipi - x'pI (89) 

up to irrelevant constants [^3| . Note that the same expression for the excess free energy of the dense system is obtained 
for a model in which the interaction energy between two species cr,a' varies as (a — a')^: 

f^^x' Jdada'p{a)p{a'){a^a'f^~x'pl+ J da p{a)a^ (90) 

The second term on the r.h.s. is linear in p(a) and can be discarded for the purposes of phase equilibrium calculations. 
The model (|8^) can therefore also be viewed as a simplified model of chemical fractionation, with a being related 
to aromaticity, for example, or another smoothly varying property of the different polymers. Eq. ( |90| ) suggests that 
the model should show fractionation into an ever-increasing number of phases as x' is increased: only the entropy of 
mixing prevents each value of a from forming a phase by itself, which would minimize this form of /. It is therefore 
an interesting test case for the moment free energy approach (and the method of adding further moment densities), 
yet simple enough for exact phase equilibrium calculations to remain feasible |^-^, allowing detailed comparisons to 
be made. 

As a concrete example, we now consider phase separation from parent distributions of the form p^^'^(a) oc exp(7(T) 
(for — 1 < cr < 1, otherwise zero). The shape parameter 7 is then a fixed function of the parental first moment density 
Pi°'' ~ ^daap^^\a). Fig ^ shows the exact coexistence curve for pf^ — 0.2, along with the predictions from our 
moment free energy with n moment densities {pi = Jdaa''p{a), i = 1. . .n) retained. Comparable results are found 

for other p^^^- Even for the minimal set of moment densities (n ~ 1) the point where phase separation first occurs 
on increasing x' is predicted correctly (this is the cloud point for the given parent). As more moments are added, 
the calculated coexistence curves approach the exact one to higher and higher precision. As expected, the precision 
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FIG. 9. Coexistence curves for a parent distribution with p\ — 0.2. Shown are the values of pi of the coexisting phases; 
horizontal lines guide the eye where new phases appear. Curves are labeled by n, the number of moment densities retained in 
the moment free energy. Predictions for n = 10 are indistinguishable from an exact calculation (in bold). 

decreases at high x', where fractionated phases prohferate; in this region, the number of coexisting phases predicted 
by the moment method increases with n. However, it is not always equal to n + 1, as one might expect from a naive 
use of Gibbs' phase rule; three-phase coexistence, for example, is first predicted for n = 4. In fact, for fractionation 
problems such as this (but not more generally), study of the low temperature limit (large x') suggests that to obtain 
n + 1 phases one may have to include up to 2n moments. We show below that using localized weight functions (rather 
than powers of a) for the extra moments can reduce this number, but also gives less accurate coexistence curves. 

In Fig. |l^, we show how, for a given value of the lever rule is satisfied more and more accurately as the number 
of moment densities n in the moment free energy is increased: the total density distribution ptot(f) — p*-"-* (cr) 

approaches the parent p^^\a). This is by design, because the moment method forces the moment densities pi = 
J da Wi{a) ptoticr) of ptotCc) to be equal to those of the parent, for i = 1 . . . n; as n is increased, therefore, ptotlc) 
approaches p'^^^{a). 

The above results for the dense copolymer model raise two general questions. First, how large does n (the number 
of moment densities retained) have to be for the moment method to give reliable results for the coexisting phases 
(and, in particular, the correct number of phases)? And second, how should the weight functions for the extra 
moment densities be chosen? Regarding the first question, note first that the theoretical framework developed so 
far only says that the exact results (and therefore the correct number of phases) will be approached as n gets large; 
it does not say for which value of n a reasonable approximation is obtained. In fact it is clear that no "universal" 
(problem-independent) value of n can exist. The model studied in this section already provides a counterexample: 
As x' is increased, the (exact) number p of coexisting phases increases without limit; therefore the minimum value 
of n (n — p — 1) required to obtain this correct number of phases from a moment free energy calculation can also be 
arbitrarily large. 

Nevertheless, we see already from Fig. |^ that, for the current example, the robustness of the results to the addition 
of extra moment densities provides a reasonable qualitative check of the accuracy of the coexistence curves. To study 
the convergence of the moment free energy results to the exact ones more quantitatively, especially in cases where 
the latter are not available for comparison, we need a measure of the deviation between ptot (tj) and p'-"H<^)- To 
obtain a dimensionless quantity that does not scale with the overall density of the parent, we consider the normalized 
distributions ntot(o') = /Otot(o')/ptot and n'-^'>{a) — p(°)(cr)/p(") and define as our measure of deviation the "log-error" 




(91) 



This becomes zero only when ntot(c) — n^^'>{a); otherwise it is positive. When the deviations between the two 
distributions are small, the logarithm becomes to leading order the relative deviation ritot(o')/n^*'-'((T) — 1 and we can 
identify S^^"^ as the root-mean-square relative deviation. We work with the logarithm rather than directly with the 
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FIG. 10. The lever rule in phase coexistences calculated from the moment free energy. We consider the same scenario as in 
Fig. ^ and take a horizontal cut through that figure, varying n, the number of moment densities retained in the moment free 
energy, at fixed x' ~ 10. (a) Solid and dashed lines: The two daughter phases calculated for n — 1 (both are exponential, but 
with different signs for the exponent). Dotted line: The corresponding total density distribution ptot{cr) = X/a=i 2 
This deviates (rather drastically, in this case) from the parent p'"' (a) (in bold) ; the moment calculation with n = 1 enforces only 
that the means pi of the two distributions be equal. (Their normalizations po are trivially equal, because we are considering 
the dense limit po = 1-) (b) Total density distributions ptot(o") for increasing values of n. Curves are labelled by the value of n; 
the curves for n = 8 and n = 10 (not labelled) are almost indistinguishable from the parent (in bold). As n increases, ptot(o") 
converges to the parent, being forced to agree with it in an increasing number of its moments, (c) Density distributions of the 
four daughter phases at x' ~ 10. For each, both the exact result and that of the moment free energy calculation with n = 10 
are shown; they cannot be distinguished on the scale of the plot. 

relative deviation because the former gives more sensible behavior for larger deviations; in particular, isolated points 
where ntot{cr) is nonzero while n^'^\a) is close to zero do not lead to divergences of S jsj]. 

Fig. |ri| shows the n-dependence of the log-error S for three fixed values of with the same parent as in Fig ^ [ p5[ . 
The correct number of phases (3, 4, 5, respectively) is reached at values of (5 « 4 x lO^'^, 5 x 10^'*, 2 x lO^'^ for x = 5, 
10, 15, respectively. This suggests that the log-error S might generally be a useful heuristic criterion for guiding the 
choice of n, with values of 6 around 10~^ (corresponding to an average deviation between ntot(o') and n*^"-* {a) of around 
1%) ensuring that the correct number of phases has been detected. Note that, although for the current example we 
did know the correct number of phases (because the simplicity of the problem made an exact calculation feasible), the 
definition of S does not require us to know in advance any properties of the exact solution of the coexistence problem. 
Also, note that in real physical materials it is extremely unlikely for the parent or "feedstock" distribution p'^^\a) 
to be known to better than 1% accuracy. Hence even for systems, such as polymers, where a mean-field (and thus 
truncatable) free energy model is thought to be highly accurate, deviations from the lever rule, at about the 1% level, 
will almost never be the main source of uncertainty in phase diagram prediction. 

We now turn to the second general question, regarding the choice of weight functions for the extra moment densities. 
Comparing ( |6l| ) with the formally exact solution ( |60| ) of the coexistence problem tells us at least in principle what 
is required: the log-ratio In i?(cr) /p^"' (ct) of the "effective prior" and the parent needs to be well approximated by a 
linear combination of the weight functions of the extra moment densities. However, the effective prior is unknown 
(otherwise we would already have the exact solution of the phase coexistence problem), and so this criterion is of little 
use Q. 

To make progress, let us try to gain some intuition from the example at hand (dense random copolymers). So far, 
we have simply taken increasing powers of a for the extra weight functions. On the other hand, looking at Fig. |l^ one 
might suspect that localized weight functions might be more suitable for capturing the fractionating of the parent into 
the various daughter phases. We therefore tried the simple "triangle" and "bell" weight functions shown in Fig. p^(a). 
For a given n, there are n — 1 extra moment densities (the first moment density, with weight function wi[a) = a, is 
prescribed by the excess free energy) and we take these to be evenly spaced across the range cr = — 1 . . . 1, with overlaps 
as indicated in Fig ^(a). (For the case of triangular weight functions, the log ratio \n p{a) / p^'^^ (a) — J^i^i'^ii'^) 
between the density distribution of any phase and that of the parent is thus approximated by a "trapezoidal" function 
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FIG. 11. Log-error 5 as a function of n for three fixed values of x' ~ 5, 10, 15 (bottom to top), for tlie same scenario as in 
Fig. 1^. The number of phases is indicated next to the curves. Note that in all three cases the correct number of phases is 
reached for 5 of order 10~^. 

across n equal intervals spanning —1 ... 1.) Fig. |l2|(b) shows the corresponding coexistence curves for the triangle 
case with n = 1,2,3,4. Note that multi-phase coexistence is generally detected for smaller n than for the case of 
power-law weight functions (Fig. |^), in agreement with the expectations outlined above. However, the predicted 
number of phases does not vary monotonically with n: in the range of x' shown in Fig. |l2|(b), up to four phases are 
predicted for n = 3, whereas n = 4 gives no more than three. This effect illustrates a significant difference between 
power-law and triangle weight functions. The former, but not the latter, make up an incremental set: increasing n 
by one corresponds to adding a new weight function while leaving the existing set of weight functions unchanged. 
The family ( [40|) of accessible density distributions p(cr) is thus enlarged without losing any of its former members. 
Intuitively, one then expects that the number of phases cannot decrease as n is increased. But for nonincremental 
sets, such as triangle weight functions (and bells, which give qualitatively similar results to triangles), each increase 
in n creates a larger family but completely discards the previous one, so that there is then no reason for the results 
to be monotonic in n. 

We regard the fact that the number of phases can go up and down with increasing n as a disadvantage of non- 
incremental weight functions. In the current example, we also find that they generally give higher values of the 
log-error than power-law weight functions; see Fig. ^(c). Finally, incremental sets of weight functions also have 
the benefit that the results for n can be used to initialize the calculation for n + 1, thus helping to speed up the 
numerics (though we have not yet exploited this property). On the basis of these observations, we would generally 
recommend the use of incremental sets of weight functions for extra moment densities. Beyond this, the optimal 
choice of extra weight functions remains largely an open problem. One avenue worth exploring in future work might 
be the possibility of choosing extra weight functions adaptively. One could imagine monitoring the log-error in a 
phase coexistence calculation, and increasing the number of extra moments by one every time it becomes larger than 
a certain threshold (lO^**, say). The new weight function could then be chosen as a best fit to the current lever-rule 
violation In ptot(f)/p^'' -'(cr) (for example as a linear combination of weight functions from some predefined "pool"). 
We are planning to investigate this method in future work. 

To conclude this section on the dense random copolymer model, we briefly discuss the spinodal criterion, and 
ask whether critical or multi-critical points can exist for a general parent distribution p'-'^)(cr) (with p^Q^ = !)• The 
criterion (^), applied to our one-moment free energy /m(pi) becomes simply 

/,„(pi) = o((pi-pl"^)'+i) (92) 

where I = 2 for a spinodal and I = 2n — 1 for an n-critical point. From eq. (|8^), we have that 

/m = Ao + Xipi - x'Pi = -Sm - x'Pi (93) 
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FIG. 12. The effect of tlie weigfit functions for extra moment densities on coexistence calculations, (a) The "triangle" and 
"bell" weight functions, for the case n — 1 = 3. (b) Coexistence curves for the triangle weight functions, for the same parent as 
in Fig. 1^. Note that multi-phase coexistence is generally detected for smaller n than in Fig.^, where power-law weight functions 
were used for the extra moments. On the other hand, the predicted number of phases now no longer varies monotonically with 
n. (c) Dependence of log-error 5 on n for fixed x' ~ 10 and three choices of weight functions for the extra moment densities: 
Power-law (solid line), triangle (dashed), bell (dot-dashed). The number of phases is marked next to the curves. 
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The excess part — xVi is quadratic in pi and therefore only enters the spinodal criterion; the additional conditions for 
critical points only depend on the moment entropy Sm = — Aq — Xipi- Because of the dense limit constraint pi = 1, 
Ao is related to Ai by 

-Ao = In / dCTp(°)(cr)e^i'" = h{Xi) 



As expected from the general discussion in Sec. [IC, — Aq is thus the cumulant generating function h{\i) of p'^^^a) 
{= n^^'' (a) in the dense limit we are considering), introduced in eq. (^) within the combinatorial approach. The 
moment entropy 3^ = h — Xipi is its Legendre transform; its behavior for pi w p^"^ reflects that of h for small Ai. We 

can thus use the expansion h = X)j!Li('^j/iOM' where the Cj are the cumulants of ^'"^(0-); ci = pf ^ is its mean, C2 
its variance etc.. Using the properties of the Legendre transform, it then follows that the leading term of the moment 
entropy is = (pi — p^'^-')^/(2c2). Inserting this into (|9^), we see immediately that the spinodal condition becomes 

2x'c2 = 2x'(pf -(pr)^)=l (94) 

To find the additional conditions for critical points, consider what happens if the cumulants C3, . . . , c; are zero, 
hence h = pJ^^Ai + C2A2/2 + Oi\[+^). It can then easily be shown that the moment entropy behaves as Sm — 

(pi — p^^^Y / (2C2) + C((pi — Pi"^)'"''"^); conversely, one can show that this behavior of Sm implies that the cumulants 
C3, . . . , c; are zero. We thus have the simple condition that for an n-critical point to be observed, the cumulants C3, 
. . ., C2„-i of the parent distribution have to be zero; the value of x' at this critical point is given by the spinodal 
condition (|9^). (For an ordinary critical point, n — 2, the condition C3 = is well known in this context [^.) We 
therefore come to the somewhat surprising conclusion that even in a very simple polydisperse system such as this - 
with a single moment density occurring in the excess free energy - multi-critical points of arbitrary order can occur, 
at least in principle. In practice, the required fine-tuning of the parent will of course make experimental study of 
these points difficult. 



2. Copolymer with solvent 

We now consider the random copolymer model in the presence of solvent, i.e., for a copolymer volume fraction 
pi < 1 . We are not aware of previous work on this model in the literature, but will briefly discuss below the link to 
models of homopolymer/copolymer mixtures ||57[| . The excess free energy ( p7| ) then depends on two moment densities, 
rather than just one as in all previous examples. For simplicity, we restrict ourselves to the case of a neutral solvent 
which does not in itself induce phase separation; this corresponds to % = 0, making the excess free energy 

/ ^ r(l - po) ln(l - po) - x'pI (95) 

This is still sufficiently simple that the exact spinodal condition ( |5^ ) - for a system with a general density distribution 
p(cr) - can be worked out analytically; one finds 



2X' = , (96) 

r[p2Po - Pl) + P2{1 - Po) 



In the dense limit po — > 1, the earlier result (p4|) is recovered. To study critical and multi-critical points, we use the 



moment free energy corresponding to (95) 



/m = Aopo + Aipi - Po -I- r(l - Po) ln(l - po) - x'pi (97) 

and start from the criterion (^4|). The moment chemical potentials are given by 

Mo = Ao - r[ln(l - po) - 1] 
pi = Ai - 2x'pi 



and the curve p(e) referred to in (54) - along which the critical phases merge - is simply parameterized in terms of 
Ao(e) and Ai(e). The restriction that the curve must pass through the parent at e = means that Ao(0) — Ai(0) — 0, 
so to get the small e behavior of the chemical potentials ( p8| ) we can expand for small Ao and Ai. To simplify the 
algebra, we restrict ourselves to parent distributions p^°^(cr) which are symmetric about cr = 0. Using this symmetry 
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FIG. 13. Coexistence curve for copolymer with solvent; the parent has a uniform density distribution 
p'°'((t) = const. = /2 with overaU density (i.e., copolymer volume fraction) Pq°' — 0.6. The results were obtained from the 
moment free energy with n = 11 moment densities retained (weight functions Wi{a) — a'', i = . . . 10). Note the tricritical 
point at x' = 2.5: The parent splits continuously into three phases. 



and the fact that e is simply a dummy parameter which can be reparameterized arbitrarily, we can then set Ai = e 
and Ao = aae^ + 046^ + O(e^). 



To have a spinodal, the moment chemical potentials (98) must differ from those of the parent by no more than 

0{e^). We find that this gives the condition 2x' = l/p'^\ in agreement with (96): symmetric parents have pf ^ = 0. 
For a critical point, the chemical potential difference must be decreased to Oie^. This gives no new condition (but 
determines 02), implying that parents with symmetric distributions are automatically critical at their spinodal. This 
result has its origin in the symmetry a — > —a of the excess free energy, and can be understood by considering the 
POj Pi: x' phase diagram of the moment free energy (|9^). In this phase diagram there is a two-dimensional spinodal 
surface, and within this surface lies a line of critical points. The symmetry a —a (corresponding to pi —pi, 
Ai —^ — Ai) forces this critical line to lie within the plane pi = of symmetrical distributions, implying that for such 
distributions spinodal and critical points coincide. 

For a tricritical point, finally, where the chemical potential difference must be further reduced to 0{e^), we find 
after some algebra the condition 



3r(pf) 



(0) 

rPo 



(0) 

Po 



As an example, consider a uniform parent distribution p^^^[a) — const, which can be written as p'^'^-'(cr) = p^^' /2, using 

JO) 



the fact that Pq"'' — J^^da p^o)(^uy Then ^2°'' = Po^V^ and pf^ = Po^V^ and so a tricritical point occurs if the overall 
density {i.e., the copolymer volume fraction) is p'^'' — 3/(2r + 3). Fig. |l^ shows the coexistence curve calculated for 

this parent (with r — 1), which clearly shows the tricritical point at the predicted value x' = 1/(2/02°'') = r + S/2 = 2.5. 
Our numerical implementation manages to locate the tricritical point and follow the three coexisting phases without 
problems; we take that as a signature of its robustness |5^. Note that the tricritical point that we found is closely 
analogous to that studied by Leibler [ST) for a symmetric blend of two homopolymers and a symmetric random 
copolymer that is, nonetheless, chemically monodisperse (in the sense that cr = for all copolymers present). In fact, 
in our notation, the scenario of Ref. simply corresponds to a parent density of the form p(°)((T) ^ S{a — 1)+S{(j+1), 
with the copolymer (cr = 0) now playing the role of the neutral solvent. 

We conclude this section with another illustration of the geometrical intuition which the moment free energy can 
provide. In Fig. 11^, we plot the moment free energy /m(po: Pi)j eq. (p7|), for the uniform parent p^°^(cr) studied above. 
Two values of x'are shown: the tricritical value (%' = 2.5), and a value just in the three-phase region (x' = 2.65); the 
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FIG. 14. Moment free energy /m(po, pi) near the tricritical point, for the same parent as in Fig. For ease of visualization, 
the free energy surface is shown upside down, with /m increasing in the downward direction, (a) At the tricritical point 
x' = 2.5, the parent is stable (note that linear terms have been added to the free energy to make the tangent plane at the 
parent horizontal), (b) As x' is increased into the three-phase region (here x' = 2.65), the parent sphts continuously into three 
phases. 



developing three minima of the free energy can clearly be seen. Even though the underlying polydisperse system has an 
infinite number of degrees of freedom p(cr), the moment method thus gives us a tool for understanding and visualizing 
the occurrence of the tricritical point in terms of a simple free energy surface with two independent variables. 



VI. CONCLUSION 



Polydisperse systems contain particles with an attribute a which takes values in a continuous range (such as 
particle size in colloids or chain length in polymers). They therefore have an infinite number of conserved densities, 
corresponding to a density distribution p{a). The fact that the free energy can depend on all details of p{a) makes the 
analysis of the phase behavior of such systems highly non-trivial. However, in many (especially mean-field) models 
the excess free energy only depends on a finite number of moment densities, i.e., (generalized) moments of p((t); the 
only dependence on other details of p{a) is through the ideal part of the free energy. For such models, which we call 
truncatable, we showed how a moment free energy can be constructed, which only depends on the moment densities 
appearing in the excess free energy. Two possible approaches for the c onstr uction of the moment free energy were 
described in Sec. ||. The fir st, w hich we call the projection method (Sec. [lA), is primarily geometrical; the second is 
mainly combinatorial (Sec. [IB) and starts fr om fi rst principles to derive an expression for the ideal part of the free 
energy in terms of moment densities. In Sec. |ll C| , we showed that these two approaches give essentially equivalent 
results. 

In Sec. [II, we explored the properties of the moment free energy in detail. The moment free energy depends only 
on a finite number of moment densities, and can be used to predict phase behavior using the conventional methods 
known from the thermodynamics of finite mixtures. We showed that, for any truncatable model, this procedure yields 
the same spinodals, critical and multi-critical points as the original free energy, even though the latter depends on 
all details of the density distribution p{a). It also correctly predicts the onset of phase coexistence, i.e., the cloud 
points and corresponding shadows. Beyond the onset of phase coexistence {i.e., when one or more coexisting phases 
occupy comparable volumes), the results from the moment free energy are approximate, but can be made arbitrarily 
accurate by retaining additional moment densities in the moment free energy. We also explained how to check the 
local and global stability of a phase coexistence calculated from the moment free energy, and how the moment free 
energy refiects overall phase behavior in the space of density distributions p{cr). 

We discussed aspects of the numerical implementation of the moment free energy method in Sec. IV, highlighting its 
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many advantages in terms of robustness, efficiency and stability. The method allows violations of the lever rule, but 
then ensures that these be kept small by adding n — K extra moments (where K is the number of moment densities 
actually present in the excess free energy of the underlying truncatable model). The solution of the phase equilibrium 
conditions thus proceeds in an n-dimensional space; but global (tangent plane) stability can always be verified by 
minimimization (of the tangent plane distance) in a space of dimension K. For K > 1 exact minimization is not 
practicable, but standard search methods can instead be employed; we chose a Monte-Carlo type algorithm for this. 

We then surveyed several sample applications in Sec. For length-polydisperse homopolymers (treated within 
Flory-Huggins theory), we showed explicitly how the moment free energy with only a single moment density (the 
volume fraction of polymer) gives the correct spinodals and critical points. By retaining the chain number density as 
an additional moment density, we could rationalize the behavior of cloud point and shadow curves as cuts through 
conventional two-species phase diagrams. As a second example, the case of copolymers with chemical polydispcrsity 
was considered, both without and with solvent. This gave us the opportunity to test the performance of the moment 
free energy method for multi-phase coexistence, and in particular to study how the accuracy of the predictions in the 
coexistence region increases as more and more moment densities are retained. In the presence of solvent, the same 
theory predicts a tricritical point under appropriate conditions, which was detected without problems, along with 
the associated three phase region, by our numerical algorithm. A plot of the two-moment free energy gave a simple 
geometrical picture of the origin of this tricritical point. 

The results for the copolymer model showed that the log-error (91) can be used to decide when enough moment 
densities have been included to locate the correct number of phases and their compositions. For general purposes we 
suggest a log-error tolerance of order lO"'^; this corresponds to limiting the rms lever-rule violations at the predicted 
phase coexistence (found by global minimization of the moment free energy) to 1% or so of the parent distribution. 
Given that in practice the parent is itself not known to such high precision, these violations lie well within an 
acceptable range of error. Indeed, every state of coexistence actually predicted by our algorithm represents the exact 
phase equilibrium of some or other parent, indistinguishable from the real one to within experimental accuracy. Here, 
as elsewhere in this paper, "exact" means exactly as predicted by a truncatable model free energy; such models are, 
of course, themselves approximate, which is another reason not to expend resources reducing the log-error below the 
level suggested above. 

In summary, with exact results for cloud point and shadow curves, critical points and spinodals, as well as refinably 
accurate coexistence curves and multi-phase regions, the moment free energy method allows rapid and accurate 
computation of the phase behavior of many polydisperse systems. Moreover, by establishing the link to a projected 
free energy fm{pi) which is a function of a finite set of conserved moment densities pi, it restores to the problem 
much of the geometrical interpretation and insight (as well as the computational methodology) associated with the 
phase behavior of finite mixtures. This contrasts with many procedures previously in use for truncatable polydisperse 
systems P-|5|,p|,p^ . Some previous approximations to the problem have used (generalized) moments as thermodynamic 
coordinates; see e.g. The moment free energy method provides a rational basis for these methods and, so 

long as the generalized moments are correctly chosen, guarantees that many properties of interest are found exactly. 
Note also that the commonplace method of "binning" the cr-distribution into discrete "pseudo-components" can be 
seen as a special case of the moment method in which each weight function is zero outside the corresponding bin; 
but since these moments do not, in general, coincide with those contained in the excess free energy the advantages 
described above are then lost. 

We finish with some remarks about extensions to the moment method, within and outside the sphere of phase 
diagram prediction. First, methods of this type may extend to nontruncatable models, for which the excess free 
energy cannot be written directly in terms of a finite number of moments as in (^. For example, many mean- field 
theories correspond to a variational minimization of the free energy: F < (E)o — TSq, where the subscript refers to a 
trial Hamiltonian | |59[ |. In such a case, one might choose to first make a physically motivated decision about which (and 
how many) moment densities pi to keep, and then include among the variational parameters the "transverse" degrees 
of freedom of the density distribution p{a). We have not pursued this approach, though it may form a promising basis 
for future developments. Secondly, we note that physical insight based on moment free energies may increasingly play 
a part in understanding kinetic problems. For example, one of us has argued that in many systems the zeroth 
moment (mean particle density) can relax more rapidly (by collective diffusion) than can the higher moments (by 
intcrdiffusion of different species) , and that this may lead to novel kinetics in certain areas of the phase diagram. And 
finally, we point to recent work |6l[| in which the moment method has been extended, via density functional theory, to 
allow study of inhomogeneous states; this opens the way to studying the effect of polydispcrsity on interfacial tensions 
and other interfacial thermodynamic properties. 
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APPENDIX A: MOMENT (GIBBS) FREE ENERGY FOR FIXED PRESSURE 



In this appendix, we explain how the construction of the moment free energy is extended to scenarios where the 
pressure (rather than the overaU system volume) is held constant. One then describes the system by the particle 
number distribution N(a) = Vp{a); this is normalized so that its integral gives the total number of particles, N = 
/ da N{a). The relevant thermodynamic potential is the Gibbs free energy G = miny F+UV, which we now construct. 
For a truncatable system, whose free energy density obeys eqs. by definition, the Helmholtz free energy can be 
written as 



F daN{a) 



In 



N{<y) 
V 



where the excess free energy F{Ni,V) = Vf{pi) depends on volume V and on the moments Ni — Vpi = 
jda N{a)wi{a) of the particle number distribution. If = 0, a Legendre transform gives the Gibbs free energy 
of an ideal mixture: 

Gid = nmi Fid + ny = T j da N{a) In ^^^^ ^ NT j da n{a) In ^(ct) + A^T In ^ 



Here n{a) — N{a)/N — p{a)/ p is the normalized particle distribution. In the general case, we split off the ideal part 

G = Gid 

and obtain for the excess Gibbs free energy 



G = minF + nV 

V 



NT 

G = minTVTln + F(N,,V) + IW 

V IW 



This is a function of N and the Ni (as well as the fixed control parameters 11, T, which we suppress in our notation). 
We can therefore write the excess Gibbs free energy per particle, G{Ni, N)/N = g{mi), as a function of the moments 
rrii = Ni/N [i = 1 . . . K) of the normalized particle distribution n{a). Altogether, the Gibbs free energy per particle 
of a truncatable system becomes 



g[n{a)] ^ G/N ^ T f da n{a) In n{a) + T In ^ + .g(m,) 



(Al) 



This has again a truncatable structure; the excess part g "inherits" its moment structure from /. Note that because 
of the normalization of the rrii, g does not depend on the density p = pQ. It is therefore normally a function of one less 
variable than /, unless / is already independent of p (as is the case in Flory-Huggins theory, eq. (^), for example). 
The chemical potentials p,{a) follow from (|Al|) as 



/i(cr) 



SN{a) 
5n{a) 



Ng[N{a)/N] 

+ g - da'n{a') 



5g 



5n{a') 



= Tlnn(o') + Wi{a) 



dg_ 

drrii 



rin-+.9-)^m.— 



(A2) 



Multiplying by nia) and integrating over a, one has g — Jda p{a)n(a) as it should be (compare eq. (Al)). Compar- 
ing (A2) with (p6[), one deduces that 



dg_ 

drrii 



and that the density p of a phase with moments rrii is given by 



Tin 



n 



E 



mip.i. 
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To construct the moment Gibbs free energy, one can now proceed as in the constant volume case: in the ideal 
contribution, we replace In n((T) — > hi[n(a) /n^^\a)], with a normalized parent distribution n^^\a), and then minimize 
g at fixed values of the rrii. The minimum occurs for distributions n{(j) from the family 



i(cr) = n'-°\a) exp ^Aq + ^ XiW^{a)^ 



(A3) 



The extra Lagrange multiplier Ao comes from the normalization condition toq = Jdan(a) = 1. Inserting into ( |Al| ), 
one obtains the moment Gibbs free energy 



i(mi) = r Ao + ^ Xtnii 



(A4) 



where the are given by 



m. 



exp Ao + ^ A, 



Wi{a) 



J da Wi{a)n^°'' {a) exp XiWi{a)) 
f dan^^Xa) exp {J^i XiWi{a)) 



The derivatives of g^ can again be identified as moment chemical potentials and obey 

, dg 



dm, 



^TX,. 



duTli 



(A5) 



The constraint mo ~ 1 does not affect this result because the partial derivative is taken with the values of all moments 
other than rrii fixed anyway. For the same reason, the second derivatives of the ideal part of g^ are (up to a factor 
T) given by the inverse of the matrix of second-order normalized moments , by analogy with ( ^ ) . (Note that the 
first row and column of the second-order moment matrix, corresponding to z = or j = 0, need to be retained; they 

can only be discarded after the inverse has been taken.) 

From (|A^), one deduces immediately that if one of a number of coexisting phases is in the family (A2), then so are 
all others. Equality of the exact chemical potentials /i(cr) in all phases is then equivalent to equality of the quantities 



jli + TXi = Hi {i 



.K), 



TXo+g-J2 



niiHi. 



(A6) 



But from (A4,A5), these are simply the moment chemical potentials = dgra/drrii and (up to an unimportant 
additive constant) the associated Legendre transform — ^^rriifj.i. Therefore, phase equilibria can be constructed 
by applying the usual tangent plane construction to the moment Gibbs free energy, in analogy with the constant 
volume case. Note, however, that this tangent plane construction now takes place in the space of the normalized 
moments rrii, rather than that of the unnormalized moment densities pi. As pointed out above, this generically 
reduces the dimension by one. For the trivial case of a monodisperse system, / is a function of p only and g does not 
depend on any density variables at all, so that the tangent plane condition becomes void as expected. 

The exactness statements in Sec. Ill can also be directly translated to the constant pressure case. The arguments 
above imply directly that the onset of phase coexistence is found exactly from the moment Gibbs free energy: all 
phases are in the family (A3), because one of them (the parent) is, and the requirement of equal chemical potentials is 
satisfied. Spinodals and (multi-) critical points are also found exactly. Arguing as in Sec. Ill A, and using the vector 
notation of (|5^), the criterion for such points is found as 



Afi = fi{n{e)) — fi{n 



(0)) 



(A7) 



where I — 2 for a spinodal and I ~ 2n ~ 1 for an n-critical point. Again, the curve n(e) only has to pass through 
phases with equal chemical potentials and can therefore be chosen to lie within the family (A3). As discussed above 
(eq. ([A^)), the chemical potential difference on the left hand side of (A7) is then zero to the required order in e if 
the same holds for the and A/i, where h ~ g^^ — rUiPi. It may look as if the constraint on h gives one more 
condition here than in the constant volume case; but in fact A/i^ — 0{e^) for all i already implies that Agm = ©(e'"'""'^) 
and hence Aft, = ©(e'). The exact condition (A7) for spinodals and critical points is therefore again equivalent to the 
same condition obtained from the moment Gibbs free energy. 
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APPENDIX B: MOMENT ENTROPY OF MIXING AND LARGE DEVIATION THEORY 



In this appendix, we discuss some interesting properties of the Legendre transform result ( |29D for the moment 
entropy of mixing, in particular its relation to large deviation theory (LDT). 

We first note that (|2^), which we derived by taking the limit a: — > of the result ( p6| ) for general x, can also be 
obtained by a more direct route. In the limit a; — > 0, the sizes of particles in the smaller system become independent 
random variables drawn from n^'^\a)] the second phase can be viewed as a reservoir to which the small phase is 
connected. One writes the moment generating function for Vim) in the small phase as a product of xN independent 
moment generating functions of 7i'^*'^((t), and evaluates the integral over V{m) by a saddle point method |36| . 

One can also view ( p9| ) as a generalisation of Cramer's large deviation theory (LDT) Recall that this theory 
treats the "wings" of a distribution like V{m) correctly, to which the central limit theorem (CLT) does not apply [|63[ . 
However, the CLT approximation in this problem has quite a beautiful interpretation, and is worth describing sepa- 
rately. It can be derived by assuming that n^^\a) is itself a Gaussian in eq. (|2^), or just written down directly: 



V{m) ~ exp 



^ (m — m^^^Y 



where w*-"-* is the variance of the parent distribution. We see that this gives a term in the free energy of the form 

Xm — TO^"^)^ 



A/ = -Tps = pT^ 



This result states that there is an entropic spring term in the free energy that penalizes deviations of the mean size 
m from the mean size in the parent distribution. The spring constant is inversely proportional to the variance in 
the parent distribution, thus if the parent is narrower, it is harder to move away from the parental mean size [ |64[ . 
However, as indicated above, the CLT is not sufficiently accurate for our purposes. For instance it gives a finite weight 
to m < even if w{a) is strictly positive. That said, the CLT may be attractive for narrow size distributions, and 
one might make a connection with the recent results of Evans et al psf (even though the latter do not seem to rely 
on a specific shape of the distribution). We have not fully explored this avenue. 



APPENDIX C: SPINODAL CRITERION FROM EXACT FREE ENERGY 



In this appendix, we apply the spinodal criterion (|5|) to the exact free energy (13^) and show that it can be expressed 
in a form identical to (p6|). This result has been given by a number of authors |ll|Jl^ , 3 , but we include it here for 
the sake of completeness. 

Choosing for p the vector p(a) (whose "components" are the values of p{a) for all a), the spinodal criterion ( |5l| ) 
applied to (BS) becomes 



da' 



- 6p{(T ) = 2^ a„ a„ ^nc^) Spj + T ■ 



5p{a)5p{a') 



dpidpj 







(CI) 



The change of \np{a) along the instability direction (which is 5p{a)/p{a)) is therefore a linear combination of the 
weight functions Wjjcr). This means that the instability direction is within the family (|4^), consistent with our general 
discussion in Sec. Ill A. One can now rewrite (CI) in the form 



Sp{a) = -(3^w,{a) p{a) 



dpidp^ 



Spj 



(C2) 



and take the moment with the fc-th weight function Wfe(cr) to get 



6pk+pJ2 



Pki 



dptdp. 



6p, = 0. 



As promised, this is identical to the spinodal condition (^^ derived from the moment free energy, with the matrix 
multiplications written out explicitly. 
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APPENDIX D: DETERMINANT FORM OF CRITICAL POINT CRITERION 



In this appendix, we give the form of the critical point criterion (^8|) that uses the spinodal determinant Y 
from ( |5^ ) [p4 |. At a critical point, the instability direction connects two neighboring points on the spinodal. The first 
order variation of the spinodal determinant Y along this direction must therefore vanish. This gives 



SY 



^dY ^ V- 5r , 



dpij 



0. 



(Dl) 



Here the Spij are the changes in the second order moment densities along the instability direction. We can express 
these as a function of the Spi via the change in the Lagrange multipliers A^: From the definition (^), a change in the 
Xi changes the first and second order moment densities according to 



dpi 
9A, 



M, 



dXk 



Pijk- 



Hence changes in second and first order moment densities are related by 



(D2) 



Spki = ^ Pku SXi = ^ pkH (M Spj 



Inserting this into (Dl) gives the determinant form of the critical point condition 

i ij k<l 



(D3) 



Of course, if in an application of the moment free energy method one succeeds in calculating Y explicitly as a function 
of the Pi alone (rather than as a function of the pi and the pij ) , then one can determine critical points simply from 
the condition 



We have written a total derivative sign here to indicate that the dependence of Y on the pij is accounted for implicitly. 
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